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FOREWORD 

This  report  summarizes  the  analysis  work  performed  on  the  hydrodynamic 
ram  effect  at  Air  Force  Flight  Dynamics  Laboratory  under  work  unit 
43630127  (Fuel  Tank  Responses  to  Projectile  Impact  and  Fluid  Shock  Pulses). 
The  work  was  performed  during  the  period  of  August  1973  through 
September  1974.  It  is  intended  that  this  study  will  form  the  basis  for 
a comprehensive  test  program  and  for  the  development  of  a refined  and 
verified  analysis  method  to  be  accomplished  in  the  next  couple  of  years. 

In  preparing  this  report,  one  purpose  has  been  in  the  foreground,  to 
create  a document  which  will  serve  as  a single  reference  for  the 
complicated  technology  irt  this  field  at  the  current  time.  It  should  be 
reasonably  complete  in  its  description  of  all  the  significant  con- 
tributions to  hydrodynamic  ram  studies.  All  of  the  major  works  as  well 
as  selected  works  of  lesser  significance  have  been  included.  All  equations 
needed  have  been  incorporated  in  the  text,  with  some  of  the  details 
of  derivation  required  for  a full  understanding  placed  in  the  appendices. 

Four  individuals  are  due  an  especial  expression  of  appreciation 
for  their  generous  contribution  of  time  while  the  original  reports  and 
papers  were  being  prepared;  they  are:  Dr.  Robert  Ball  : of  the  Naval 

Postgraduate  School,  Dr.  Eric  Lundstrom  of  NWC,  Dr.  K.S.  Nagaraja  of  the 

AFFDL  and  Dr.  Peter  Torvik  of  the  Air  Force  Institute  of  Technology. 

Others,  whose  patience  and  encouragement  were  most  helpful  and  appreciated 

were  Mr.  Andre  Holten  and  Lt  Col  Duane  Baker  of  AFFDL.  To  all,  the 
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SECTION  I 
INTRODUCTION 

Vietnam  combat  experience  has  demonstrated  that  aircraft  can  be 
particularly  vulnerable  to  small  arms  and  automatic  fire.  Typically, 
a small  arms  round  can  penetrate  the  skin  and  fuel  tank  of  a low  fly- 
ing aircraft.  If  the  fuel  which  escapes  from  this  impacted  area  ig- 
nites, then  an  explosion  and  fire  can  bring  the  aircraft  down.  In  ad 
dition,  small  arms  rounds  can  cause  excessive  structural  damage  as 
their  kinetic  energy  is  given  up  to  the  fluid  and  then  transferred 
to  the  fuel  tank  walls.  In  the  past,  this  phenomenon  has  been  called 
hydraulic  ram,  but  this  term  already  refers  to  an  engineering  device. 
It  is  suggested  that  the  term,  hydrodynamic  ram,  is  both  more  descrip 
tive  and  more  appropriate. 

A similar  rupture  phenomenon  had  been  a problem  for  NASA  when 
hypervelocity  particles  impacted  the  liquid  oxygen  fuel  tanks  of 
spacecraft.  NASA  made  a series  of  studies  to  measure  and  explain 
the  event,  but  was  able  to  eliminate  their  problem  without  complet- 
ing the  research.  Thin  metal  bumpers,  placed  in  front  of  the  fuel 
tank,  break  the  meteorite  particles  into  a shower  of  very  small, 
light  pieces.  The  spacecraft  fuel  link  walls  can  then  withstand 
the  distributed  load. 

A similarly  simple  solution  is  not  possible  to  protect  aircraft 
against  ordnance  projectiles;  further,  it  is  not  reasonable  to  design 
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the  structure  to  withstand  the  Impact,  because  of  the  weight  penalty. 
Thus,  it  became  necessary,  from  a military  standpoint,  to  study  hy- 
drodynamic ram  more  closely  in  order  to  find  clues  to  an  efficient 
defense  mechanism.  As  a result,  several  additional  studies  have  been 
completed;  however,  in  most  cases  each  study  treats  only  one  or  two 
aspects  of  the  hydrodynamic  ram  event  in  detail.  Thus,  a complete 
grasp  of  hydrodynamic  ram  can  only  be  gained  by  studying  more  sources 
and  there  is  the  real  possibility  that  some  worthy  discussions  or  an- 
alyses might  be  omitted  simply  because  they  were  published  several 
years  ago.  A repetition  of  the  known  analyses  and  explanation  of 
hydrodynamic  ram  in  consistent  notation  is,  therefore,  necessary  for 
both  completeness  and  perspective,  and  is  the  subject  of  this  report. 

1.  THE  HYDRODYNAMIC  RAM  EVENT 

Hydrodynamic  ram  is  initiated  by  the  impact  of  a projectile  into 
a liquid-filled  container.  Detailed  analyses  explain  the  event  by 
using  different  models  along  the  trajectory.  But  it  is  obvious, 
even  without  knowledge  of  those  models,  that  the  projectile  is  going 
to  give  up  some  of  its  energy  along  every  increment  of  its  path. 

This  lost  projectile  energy  is  transferred  to  the  fluid  as  kinetic 
and  thermal  energy;  the  fluid  then  transfers  its  kinetic  energy  to 
the  tank  walls.  Thus,  hydrodynamic  ram  is  fluid  loads,  especially 
impulsive  loads,  on  the  fuel  tank  walls. 
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In  sequential  order,  the  phases  of  hydrodynamic  ram  are  impact, 
fluid  shock,  fluid  drag,  and  exit.  These  phases  are  illustrated  in 
Figure  1 . 
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The  impact  phase  occurs  as  the  projectile  penetrates  the  tank 
wall,  and  subsequently,  shock  waves  are  induced  in  both  the  wall  and 
the  projectile.  This  phase  is  very  brief,  typically  less  than  0.1 
msec.  (A  projectile,  60  irni  in  length,  traveling  at  an  average  ve- 
locity of  0.85  km/sec,  would  be  in  contact  with  the  entrance  wall  for 
51  usee.)  Because  this  phase  is  so  brief  and  difficult  to  observe, 
there  is  no  analysis  available  which  treats  specifically  the  impact 
portion  of  the  hydrodynamic  ram  problem.  It  should  be  presumed, 
however,  that  the  stresses  created  by  the  projectile  at  impact  will 
cause  small  cracks  to  form  and  will  radically  alter  the  load-bearing 
capacity  at  the  impact  region  so  that  the  ensuing  high  pressure  will 
find  the  impact  region  the  weakest  part  of  the  tank  walls.  Never- 
theless, for  simplicity  of  analysis,  it  is  normally  assumed  that  the 
projectile  will  simply  create  a hole  slightly  larger  than  the  projec- 
tile impact  area  in  the  front  wall,  and  all  other  impact  effects 
are  ignored. 

The  second  phase  (fluid  shock  phase)  of  the  hydrodynamic  ram 
event  is  characterized  by  a local,  nearly  hemispherical,  shock  sur- 
face in  the  fluid  which  is  formed  about  the  impact  point.  This  second 
phase  is  also  of  extremely  short  duration,  beginning  at  impact,  and 
lasting  less  than  0.2  msec.  Usually  the  shock  radius  reaches  a maximum 
of  about  20  cm.  It  is  easily  observed  with  shadowgraph  photography, 
and  consequently  there  is  considerable  analysis  available.  The  shock 
phase  is  reported  ir.  detail  by  NASA  researchers,  and  it  is  the  only 
phase  they  reported  since  the  light  particles  they  studied  spend 
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nearly  all  of  their  energy  on  impact.  The  significance  of  this  phase 
was  established  by  Stepka  and  Morse  (1:13)*  who  observed  that,  for 
small  projectiles  fired  through  prepunched  holes,  fracture  took  place 
between  27  and  40  msec  after  impact.  The  model  usually  replaces  the 
projectile  with  an  impact  impulse  over  a very  brief  increment  of 
time.  The  analysis  involves  determining  the  location  and  velocity  of 
the  shock  surface,  and  evaluating  the  very  high  pressures  behind  the 
shock  surface.  This  initial  severe  loading  of  the  impacted  wall  by 
the  fluid  within  the  shock  region  can  be  expected  to  produce  or  ex- 
tend fractures  in  every  case,  and  hence  to  further  modify  the  material 
properties  of  the  wall. 

The  third  phase  lasts  for  several  milliseconds;  it  is  thus  rela- 
tively long,  and  has  been  very  closely  observed  and  analyzed.  Although 
it  is  sometimes  called  the  cavity  phase,  the  term  "drag  phase”  is  pre- 
ferred, because  the  pressure  field  can  be  accurately  related  to  the 
fluid  drag  of  the  tumbling  projectile.  A gaseous  cavity  is  an  easily 
observed  characteristic  of  this  phase;  it  is  formed  as  the  projectile a 
imparts  an  outward  radial  velocity  to  the  fluid.  The  cavity  surface 
continues  to  expand,  even  after  the  projectile  has  exited  the  tank, 
its  shape  finally  becoming  nearly^spherical  (Figure  2).  The  cavity 
collapses  to  some  minimum  volume,  then  rebounds  on  itself.  This 
cycle  can  be  repeated  three  or  four  times  before  the  cavity  dies. 


* Reference  1,  page  13. 
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Though  the  pressure  field  is  not  so  easily  observed,  it  Is  the 
most  Important  characteristic  of  the  drag  phase.  Generally,  a pres- 
sure transducer  located  where  it  sees  no  reflections  will  detect  a 
single  pressure  pulse  (2),  and  this  pulse  may  have  a steep  enough 
front  that,  for  some  applications,  it  could  be  termed  a shock.  Most 
pulses  are  less  than  0.5  msec  and  some  are  spikes  of  less  than  0.1 
msec.  The  pressure  field  is  not  uniform;  much  higher  pressures 
are  observed  near  the  projectile  when  it  is  in  a fully  tumbled  posi- 
tion. Typical  peak  pressures  at  about  15  cm  from  the  bullet  are  less 
than  70  and  are  seldom  over  140  bars  from  projectiles  up  to  14.7  mm 
API.  Two  pressure  pulses  are  shown  in  Figure  3,  where  the  second 
curve  would  be  considered  typical. 
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(t=0  is  adjusted  to  beginning  of  pulse.) 
Figure  3.  Pressure  History  Plots 
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Because  the  exit  phase,  like  the  impact  phase,  is  also  very 
brief,  its  effects  have  been  neglected  in  testing  programs  and  by 
analysis.  It  is  clear  from  high  speed  photographs  that  the  projec- 
tile will  tumble  almost  invariably  in  the  fluid,  and  will  thus  exit 
in  a tumbled  condition  most  of  the  time.  Further,  since  the  tumbled 
projectile  is  preceded  by  a very  high  pressure  region,  it  is  supposed 
that  the  bullet  strikes  an  exit  wall  that  is  already  under  some 
stress.  To  support  this,  researchers  have  occasionally  noted  that 
the  rear  wall  can  be  seen  to  bulge  before  the  projectile  penetrates  it. 
However,  close  examination  and  measurement  of  this  event  are  totally 
lacking. 

2.  OBJECTIVES 

There  are  two  major  objectives  to  this  study  of  the  hydrodynamic 
ram.  The  first  is  to  find  means  of  attenuation,  the  reduction  of  the 
damage  to  some  minimum  and  acceptable  level.  But  in  spite  of  the  great 
amounts  of  data  and  complicated  models  and  analysis,  no  completely 
satisfactory  solution  has  been  forthcoming.  An  adequate  solution  to  the 
hydrodynamic  ram  rupture  is  expected  to  be  some  combination  of  materials 
and  geometry  which  will  spread  the  pressure  pulse  out  over  space  or  time 
or  both.  This  dispersion  of  the  pulse  must  be  to  such  an  extent  that 
each  element  of  the  wall  experiences  a sufficiently  small  impulse  so 
that  its  yield  limits  are  not  exceeded.  Of  the  many  materials  that  have 
been  suggested  to  reduce  the  damage,  only  foams  (flexible  and  rigid) 
have  shown  any  promise.  Although  some  benefit  from  attenuation  mater- 
ials is  anticipated,  it  appears  now  that  this  benefit  will  be  quite 
limited.  Ram  attenuation  is  discussed  in  the  last  section. 
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Since  ram  attenuation  has  proved  so  elusive,  greater  importance  has 
been  placed  on  the  other  major  objective,  predicting  the  location  and 
extent  of  the  damage.  This  objective  can  be  met  by  developing  and 
verifying  the  analytical  equations,  many  of  which  are  the  substance  of 
the  following  pages. 
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SECTION  II 

HYDRODYNAMIC  EQUATIONS 

In  this  section,  some  of  the  equations  basic  to  the  studies  of 
hydrodynamics  are  presented.  Since  their  inclusion  is  intended  only  to 
refresh  the  reader  who  already  has  an  understanding  of  hydrodynamics  and 
to  serve  as  a point  of  reference  for  later  sections,  the  presentation 
contains  very  little  explanation.  The  following  topics  are  covered: 
conservation  equations,  jump  conditions,  thermodynamic  relations,  ve- 
locity potential,  and  Bernoulli's  equation.  Some  of  these  equations  are 
then  repeated  in  nondimensional  form  following  the  work  of  Bach  and  Lee. 

1.  CONSERVATION  EQUATIONS 

The  conservation  equations  can  be  derived  by  considering  an  arbi- 
trary, infinitesimally  small  volume  which  is  allowed  to  move  with  the 
fluid.  Then  one  assumes  there  are  no  viscous  stresses,  that  heat  con- 
duction is  negligible,  and  that  there  are  no  discontinuities  of  pres- 
sure, velocity,  or  internal  energy.  Under  these  conditions,  the  mass, 
momentum,  and  energy  of  each  elemental  volume  are  conserved.  Conservation 
of  Mass  (Continuity  Equation):  The  change  in  density  is  related  to  a 
corresponding  change  in  the  dimensions  of  the  small  volume 

If  + ^v-r-o  o) 


where  p is  the  density  and  y is  the  particle  velocity. 
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M u 

Conservation  of  Momentum:  Newton's  second  law,F»m  , is  applied  to 

d T 

the  elemental  volume: 


+ _L  7p  s o (2) 

dt  p 

i 

i 

i 

i 

where  P is  the  pressure.  / 

/ 

Conservation  of  Energy:  The  work  done  on  a small  volume  of  constant 

mass  is  equal  to  the  change  in  internal  and  kinetic  energy  of  the  volume 

/ 

/>  *-[• + -¥-]  = -*  • «*»*  (3) 

where  e is  the  specific  internal  energy  (energy  per  unit  volume).  If 
now  the  mass  and  momentum  equations  are  combined  with  the  energy  equa- 
tion, the  result  is  a more  useful  form  for  the  conservation  of  energy 
equation: 


In  all  these  equations,  it  must  be  noted  that 


+ U • V) 
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since  the  volume  has  been  allowed  to  move  with  the  fluid.  One  can 
convert  from  moving  to  fixed  elemental  volume  by  using  Equation  5. 


For  shock  wave  work,  the  conservation  equations  ars  most  frequently 
written  in  spherical  coordinates.  Dy  referring  to  Appendix  A,  these 
equations  can  be  easily  written  out,  and  are  listed  in  Lagrangian  coor- 
dinates as  follows: 


Conservation  of  Mass 


i£_  + . J!a.  + u 


+ 2^.0  (6) 


du  , du 

Conservation  of  Momentum  p *T"  + P.u  ~ — 

ot  ot 


+ 4^  = 0 


Conservation  of  Enerc 


I K.  * ' 

p"  t *•-*. 

f k> 


(*  * -Sr) 


Equations  6,  7,  and  8 are  valid  throughout  the  fluid,  except  across 
a shock  front.  The  abrupt  changes  across  the  shock  surface  are  treated 
as  discontinuities,  and  the  proper  equations,  usually  referred  to  as 
jump  conditions,  must  be  found.  These  relations  can  be  derived  by 
referring  to  Figure  4 and  by  considering  what  an  observer  situated  at 
the  shock  surface  would  see  (3:33-36).  The  following  equations  then 
express  conservation  of  mass,  momentum,  and  energy  across  the  shock 
front. 
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Figure  4.  Jump  Conditions 


*0 » po ,U0 


Mass 


/>0(u«‘  uo)  = P\  (U*  “ Ul) 


7>0(u,'  »o)(ui  ~ uo)=  pi~po 


Energy 


- po«o=  p0( US 


where  subscripts  o represent  initial  conditions  ahead  of  the  shock; 
subscripts  I represent  shock  conditions  immediately  behind  the  shock 
surface;  and  U#  is  the  velocity  of  the  shock  surface. 
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Since  U0  is  nearly  always  zero,  these  equations  are  easily  simpli- 
fied for  most  applications: 


Mass  pQ U,  = (U,  - u,) 


(12) 


Momentum  p0  Uf  u,  s Pt  - 


(13) 


Energy 


(14) 


Finally,  velocity  terms  can  be  eliminated  from  Equations  12-14  to  get 
the  Rankine-Hugoniot  relation, 


(15) 


2.  THERMODYNAMIC  RELATIONS 

Several  relations  from  thermodynamics  play  a key  role  in  the  anal- 
ysis of  hydrodynamic  shocks  (3  and  4).  The  first  relation  has  already 
been  implied  in  the  assumptions  required  for  the  conservation  equations, 
i.e.,  heat  conduction  and  viscosity  can  be  neglected.  The  assumption 
that  there  is  no  heat  exchange  between  neighboring  volumes  is  the  same 
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as  assuming  that  the  changes  in  the  states  of  the  infinitesimal  volume 
are  adiabatic,  or  that  entropy  is  constant  throughout.  Constant  entropy 
throughout  the  region  is  given  by 

42.  = o (16) 

d» 

where  S is  the  specific  entropy.  A second  important  relation  states 
that  the  change  in  specific  internal  energy  is  given  by  the  change  in 
the  heat  minus  the  work  done  on  the  volume: 

de  = TdS  - Pdr  (17) 


where  T is  the  temperature  and  t is  the  specific  volume  (volume  per  unit 
mass).  Note  that  px  - I . 

Finally,  since  the  conservation  equations  give  three  equations  for 
four  unknowns,  a fourth  equation  is  needed.  An  equation  of  state  can  be 
chosen  that  relates  any  three  of  the  variables, P , T ,p  ,r  , and®  , 
because  it  is  known  from  thermodynamics  that  only  two  of  these  are 
really  independent.  For  example,  for  a perfect  gas,  the  temperature  can 
be  written  as  a function  of  the  pressure  and  volume:  PV=  RT.  We  will 
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need  the  Tait  equation  of  state  which  can  be  written  in  the  form 
P--Fip,S)  , where  F represents  the  equation  of  state.  Specifically,  we 
will  use  the  adiabatic  form  of  the  Tait  equation  of  state  PsF(^)  , given 
here  without  derivation  (3:44): 


P = 


(18) 


where  C0  = sound  velocity  in  the  undisturbed  fluid  and  n is  taken  to  be 
constant. 

3.  VELOCITY  POTENTIAL 

A number  of  simplifications  can  be  made  for  fluid  flow  which  is 
radial.  For  this  condition,  the  curl  of  the  velocity  (VXu)  is  zero, 
and  hence  the  velocity  can  be  written  as  the  gradient  of  the  scalar 
velocity  potential :u  = V<f>  . If,  in  addition,  the  fluid  velocity  is 
small  and  nearly  constant,  the  potential  satisfies  the  wave  equation. 

It  is  useful  to  define  some  additional  terms  before  showing  this 
last  statement: 

Specific  Enthalpy:  I = • + t P (19) 


Kinetic  Enthalpy: 


I * i " »0  + 


2 


(20) 
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Sound  Velocity:  C2  = ( dP  A (21) 

\ *P  /, 

where  i0  is  the  enthalpy  of  the  undisturbed  fluid  and  the  derivative  in 
Equation  21  is  taken  at  constant  entropy  (4:84). 

Finally,  by  combining  Equation  17  with  the  differential  form  of 
Equation  19,  we  obtain 


di  = TdS  + r dP 


(22) 


The  conservation  of  mass  equation  (Equation  1)  can  be  rewritten  so 
that  it  illustrates  for  irrotational  flow  (VXu*0)  the  conditions  for 
which  the  scalar  potential  satisfies  the  wave  equation.  Equation  1 
can  be  rewritten  as 


V • if 


JL  -It 

p dt 


(23) 


It  is  shown  in  Appendix  B that  with  the  assumption  of  constant 
entropy.  Equation  23  can  be  written  as 


V • if  s 


(u  • u“)  - ^ ( 7f»V)  (if  • "IT ) (24) 

<7t 
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Since  u=V^>  , for  irrotational  motion,  and  I a-  (see,  also, 

Pt 

Appendix  B),  the  assumption  of  irrotational  motion,  constant  entropy, 
and  small  velocities  and  velocity  derivatives  implies  that  ^ satisfies 
the  wave  equation 

v‘*  (26) 

4.  BERNOULLI'S  EQUATION 

Bernoulli's  equation  is  the  basic  relation  to  be  solved  to  find  a 
description  of  the  pressure  field  in  the  tank.  The  general  form  of 
Bernoulli's  equation  can  be  derived  (3:286)  by  considering  the  con- 
servation of  momentum  equation  (Equation  2)  when  the  system  is  exposed 

-A 

to  some  external  force,  F : 

~ + -L  VP  * ? (26) 

<iT  p 


The  form  of  Equation  26  is  converted  to  a fixed  elemental  volume  by 
using  Equation  5.  At  the  same  time,  only  irrotational  flow  and  con- 
servative external  forces  are  considered  (F*7il  where £2  is  a scalar  force 
potential);  thus 


di  +(«  V)TT 


+ -L  7p  : V£1 
P 
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Since  usV^.Iu’VJTT  can  be  replaced  with-^Vu*.  Interchanging  space 
and  time  derivatives,  and  assuming /o  constant.  Equation  27  can  be  re- 
written and  integrated  to  yield  the  generalized  form  of  Bernoulli's 
equation 


dt 


SI  + 1 (t) 


(28) 


where  f (t)  is  a constant  of  integration  and  hence  is  a function  of  time 
only.  Equation  28  has  been  derived  for  irrotational , incompressible 
flow  involving  only  conservative  forces.  Allowing  f(t)  to  be  zero  would 
further  require  the  equation  be  limited  to  steady  flow.  Finally,  it 
should  be  noted  that  the  terms  have  units  of  energy  per  unit  mass. 

5.  NONDIMENSIONAL  FORM 

One  of  the  tools  used  in  the  analysis  of  shock  waves  formed  by  the 
point  release  of  energy  is  the  principle  of  similarity.  This  principle 
states  that  the  properties  of  shock  waves  do  not  change,  even  though  the 
energy  released  does  change,  provided  that  the  time  and  distances  by 
which  the  event  is  measured  are  scaled  by  the  amount  of  energy  difference. 

The  advantages  of  similitude  are  that  it  allows  the  correlation  of 
results  from  different  impact  conditions,  and  that  it  permits  reasonable 
assumptions  about  the  relative  importance  of  terms,  often  permitting 
simplifications  that  might  not  normally  be  apparent. 
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There  are  two  forms  of  similarity,  the  general  form  and  the  self- 
similar form,  the  latter  being  the  simpler.  The  goal  is  to  express  the 
hydrodynamic  equations  in  terms  of  scaled  (nondimensionalized)  functions 
that  can  be  made  to  satisfy  conditions  for  a self-similar  solution. 

Thus,  scaling  factors  must  be  found  so  that,  with  varying  projectile 
impacts,  the  density,  pressure,  and  particle  velocity  profiles  will  have 
the  same  form.  The  primary  reference  for  this  section  is  Bach  and  Lee 
(5). 

Letting  the  subscript  s0  Indicate  quantities  along  the  axis  of 
symmetry  (horizontal  axis),  e0  be  the  sound  velocity  in  the  undisturbed 
region,  M#q  be  the  shock  mach  number  along  the  axis  of  symmetry,  and  the 
bar  above  the  symbol  indicate  that  it  is  a nondimensional  quantity,  the 
following  definitions  are  used  to  write  the  hydrodynamic  equations  in 
nondimensional  form. 


Shock  Mach  Number: 


Position: 


(29) 


(30) 


5 

■j 
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(31) 


(32) 
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Pressure: 


* (f. «..) 


Specific  Internal  Energy: 


Pir, t ) 

/>  U, 

o 


t |r,t) 


* (7  * M»o)  * 


(t) 


(33) 


(34) 


For  all  this  work,  the  shock  wave  is  considered  hemispherical,  and 
therefore  uSo=U,  , and  M,o=Ms  . 

The  conservation  equations  can  now  be  given  in  spherical  coor- 
dinates and  in  terms  of  the  nondimensional  variables  just  defined. 

Conservation  of  Mass: 


(»--)  4-  * 7 -4-  + 2 ^ 

df  df  r 


= - /3m, 


dp_ 

dM, 


(35) 


Conservation  of  Momentum: 


i *.  s 

H- 

i r 

'■  Vi 

‘ * 'f 


(u -T) 

Conservation  of  Energy: 


4?-  + + 4r 

dr  p dr 


= -/3m, 


(u  -T) 


h P_ 

dr  dr  J 


+ z/3s  * ~ /3m 
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du 

dM, 


(36) 


[di  p d/o 

dM,"  pi  dti. 


(37) 
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where 


U? 


(38) 


The  derivation  of  Equations  42-44  is  given  in  Appendix  C. 


The  simplified  jump  conditions  (Equations  12,  13,  and  15)  can  also 

be  written  in  nondimensional  form.  The  subscripts  I and  o mean  the  same 

as  before  so  that  a a £i-  , etc.  Also,  Equations  13  and  15  are  further 
r 

simplified  by  recognizing  that/>0  and  e0  are  small  enough  that 
p'-p0ap,  and  e,±e0ai8|  . After  these  small  simplifications,  the 

nondimensional  jump  conditions  become 


Mass: 

1 

I 

Cl 

ft 

(39) 

Momentum: 

P,  « U, 

(40) 

Energy: 

T * M 1 ' 

*> 

(41) 

Finally,  the  Tait  equation  of  state  can  also  be  written  in  non- 
dimensional  form: 


(42) 
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6.  SHOCK  FRONT  FORMATION 

A point  source  disturbance  in  an  infinite  fluid  will  create  a 
pressure  pulse  that  moves  radially  away  from  the  point  of  the  distur- 
bance and  which  decays  in  amplitude.  The  fluid  characteristics  within 
the  pulsed  region  are  not  constant;  the  second  velocity  within  the 
pulsed  region  varies  depending  on  the  pressure  and  the  density.  Since 
the  sound  velocity  is  greatest  in  regions  of  highest  pressure,  the 
disturbance  will  propagate  more  quickly  through  those  regions.  The  net 
result  is  as  shown  in  Figure  5;  the  high  pressure  disturbance  tends  to 
overtake  the  leading  lower  pressure  disturbances,  and  thus  a shock  tends 
to  form  at  the  front  of  every  pressure  pulse.  Ana,  continuing  this 
argument,  the  back  part  of  the  pulse  will  tend  to  spread  out. 


Figure  5.  Shock  Front  Formation 
23 
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SECTION  III 
FLUID  SHOCKS 

A shock  front  is  a region  of  the  fluid  characterized  by  abrupt 
changes  in  fluid  variables  such  as  pressure,  density,  and  particle 
velocities.  Fluid  shocks  are  invariably  created  at  impact  and  are 
frequently  observed  at  points  adjacent  to  the  projectile  trajectory. 

The  latter  shocks  dissipate  quickly  as  the  disturbed  region  expands 
radially  and  thus  have  not  been  carefully  observed  or  analyzed.  Fluid 
shocks  at  the  impact  point,  however,  have  been  extensively  photographed, 
measured,  and  analyzed,  and  are  the  subject  of  this  section. 

The  impact  shock,  formed  as  the  projectile  strikes  the  tank  wall, 
creates  a hemispherical  shock  surface  which  expands  about  the  impact 
point.  If  the  projectile  is  a very  light  particle,  it  may  not  pene- 
trate the  tank  wall,  or  it  may  just  penetrate  and  be  essentially 
stopped  within  a few  centimeters.  A heavier  projectile,  such  as  a 
bullet,  can  be  expected  to  penetrate  the  wall  and  have  sufficient 
velocity  remaining  to  continue  through  the  fluid,  transferring  some 
of  its  kinetic  energy  to  the  fluid  at  every  increment  along  its  path. 

In  the  first  case,  the  impact  shock  is  due  to  a point  source  release 
of  energy,  and  it  is  this  phenomenon  that  was  observed  and  reported 
by  NASA.  In  the  second  case,  the  projectile  cannot  be  a point  source 
of  energy  except  in  the  limiting  case  or  as  an  approximation.  Since 
accurate  measurements  of  the  variables  associated  with  impact  shock 
of  a bullet  have  not  been  made,  this  section  reports  observations 
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made  by  NASA  of  light  hypervelocity  particles.  This  is  followed  by  the 
shock  phase  analysis  of  Yurkovich  which  was  intended  to  be  used  for 
ordnance  projectiles,  but  has  only  been  verified  against  data  from  light 
hypervelocity  particles. 

1.  DESCRIPTION  OF  THE  FLUID  SHOCK  ABOUT  THE  IMPACT  POINT 

In  their  preliminary  investigations,  Stepka  and  Morse  (1)  studied 
front  wall  fracture.  They  fired  small  spheres  and  cylinders  (nylon, 
aluminum,  steel,  copper,  and  tungsten)  into  aluminum  tank  wall  specimens 
mounted  on  a transparent  plastic  tank.  Some  of  the  projectiles  were 
fired  into  specimens  with  pre-punched  holes  to  eliminate  the  effects  of 
the  wall -projectile  interaction.  The  shock  waves  were  monitored  by 
piezoelectric  crystal  pressure  sensors  near  the  impact  point  (36  6 mm 
and  47.5  mm)  and  by  a high-speed  continuous-writing  camera  throughout 
the  fluid  (up  to  1.6  frames/ysec) . These  were  comparatively  low-energy 
experiments,  the  kinetic  energy  of  most  of  the  shots  being  between  100  J 
and  1350  J.  The  results  most  significant  to  hydrodynamic  ram  investi- 
gations are  that  velocities  were  clearly  an  important  factor  in  causing 
fracture,  and  that  tank  volume,  liquid  level,  and  compressibility  were 
not.  Further,  it  was  found  that  pressurization  of  the  tank  reduced 
fracturing  for  high  velocities.  It  was  also  noted  that  fracture  oc- 
curred between  27  and  40  psec  after  impact,  while  the  shock  front  was 
still  only  several  centimeters  in  radius.  Pressures  of  approximately  7 
kilobars  were  calculated  at  about  15  mm  from  the  impact  point.  Figure  6 
shows  the  pressure  sensor  records.  Finally,  it  can  be  seen  from  Figure 
7 that  fracture  thresholds  exist;  for  a given  projectile  and  target,  the 
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A.  Pressure  measured  36.6  mm  from  impact  point 

B.  Pressure  measured  47.5  mm  from  impact  point 

Impacting  particle:  aluminum  sphere,  5.55  mm  diameter; 

velocity  - 2.039  mm/ysec 

Figure  6.  Pressure  Histories  Near  the  Region  of  Hypervelocity  Impact 
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Target:  1.59  mm  aluminum  wall  of  a water- 

Projectile:  spheres;  Mass(gm) 

-filled  tank 
Diameter (mm) 

Material 

Symbol 

1.34 

5.55 

Tungsten 

A 

.699 

5.55 

Steel 

o 

.251 

5.55 

Aluminum 

□ 

.105 

5.55 

Nylon 

o 

.249 

3.18 

Tungsten 

bk 

.130 

3.13 

Steel 

V 

.019 

2.33 

Aluminum 

A 

.016 

1.59 

Steel 

< 

Open  symbols:  puncture  only 

Solid  symbols:  fracture/rupture 
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type  of  damage  changes  from  puncture  to  fracture  as  the  velocity  of  the 
projectile  is  increased. 

In  a follow-on  report,  Stepka,  Morse,  and  Dengler  (6)  investigate 
the  characteristics  of  the  pressure  waves  in  liquid-filled  tanks.  A 
test  set-up  similar  to  that  in  the  preliminary  investigations  was  used 
with  the  addition  of  a Kerr  cell  shadowgraph  system  and  a new  gun  capa- 
ble of  velocities  in  the  2.7  mn/ysec  to  6.4  mm/ysec  range.  The  primary 
goal  was  to  measure  the  shape,  velocity,  and  the  time-rate-of-change-of- 
magnitude  of  the  shock  wave. 

The  shape  of  the  wave  was  determined  to  be  hemispherical  regardless 
of  the  obliquity  of  the  impacting  particle.  McMillan  had  reported 
earlier  (7:205)  that  the  shock  surface  was  actually  ellipsoidal  with 
very  small  eccentricity,  and  was  with  the  longer  axis,  the  one  in  the 
direction  of  the  projectile  trajectory.  Velocities  of  the  shock  front 
can  be  determined  frcm  plots  of  the  shock  position  as  a function  of  time 
(Figure  8)  and  it  can  be  seen  that  the  slopes  of  these  curves  can  be 
approximated  by  straight  lines  on  a log-log  plot.  Thus 

A log  r 

Alogt  ' a (43) 


or 


(44) 


Srti’flEfer  -j 
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where  a is  a nearly  constant  slope.  The  exponent  a is  between  0.7  and 
0.9;  for  the  copper  cylinder  in  Figure  8, a is  0.733  between  10  and  20 
ysec  after  impact. 

The  projectile  deceleration  is  indicated  in  Figure  9.  The  same 
relationships  hold  as  before,  but  with  the  slope  a not  nearly  so 
constant.  However,  for  the  more  massive  particles,  as  in  tests  10  and 
11,  a is  basically  constant  at  about  0.77  and  0.91,  respectively. 

The  maximum  pressures  occur  just  behind  the  shock  surface,  immed- 
iately after  impact,  and  are  shown  in  Figure  10.  These  measurements  are 
not  made  directly,  but  shadowgraph  photography  is  used  to  learn  the 
velocity  of  the  shock  front.  Shock  front  velocity  is  then  related  to 
pressure  at  the  shock  front  through  the  Tait  adiabatic  equation  of 
state,  the  specific  volume,  and  the  specific  heat,  as  plotted  in  Figure 
11.  Note  from  Figure  10  that  the  maximum  pressures  behind  the  shock 
front  decay  vary  rapidly  with  distance.  Thus,  for  the  highest  impact 
energies  shown  here  (8254  J),  and  presumably  for  considerably  higher 
impact  energies,  the  extremely  high  pressure  region  remains  local. 

Thus,  damage  to  the  side  and  rear  walls  is  not  normally  related  to  the 
fluid  shock  produced  at  front  wall  impact. 

The  energy  of  the  projectile  can  be  measured  at  impact,  and  if 
clear  photographs  are  available,  its  energy  after  impact  can  be 


AFFDL-TR-75-102 


lO  O H ^ fO 


h rs  tn 
n co  ro 


UIQ  CtL  ^ 
c*:  »—.  uj  ct 

UJ  CQ  O UJ 

ia  z x uj 

n_  «q:  h o.  a 

00  O _J  oo  w 
>-  _ X 

s:  z o ac  o. 
X Sy  = U1 
2:  i-c;  o;  z 

_ £,UJ  UJ  ►-.  _j 
£ 03  3:  O.  SE  UJ 
32CL  fl.  X UJ 

_j  g</>  OJH 
<h  O <C  00 


C O 

a;  +->  aj 

E C >- 

0»*-J  O 
(O  3: 

i-  a> 

Ll  E'o 

•r-  dJ 

<u  t-  jz 

<—  OJ 
■r-  (/>  c 

+j  > 3 

O CL 
0)  -C  OJ 

•o+*  i- 

<_>  O QlO- 

uj  s-  oj 

a)  a.  o x: 

— s — <•  cn 

^ c 3 

A 0)0  0 

LU  r-  *r-  J- 

2:  ••-+->  -C 

M -P  <D  H 

| — <->  i- 

(U  P i. 

•r^  0)  a) 

o c -P> 

i-  a;  fu 

D.  Q.  3 


in 

VO  lO  i£)  N 

ID  If)  ift  rt  id 

mm  m co  f-t 

OD  O > v 


000 

o co 


(WW)  3DNV1SIQ 


H-  h-  O 
UwUI 
<UID 
CLO'N 


tnr^  (vu 

CO  C\J  CO  < 


cc  a o.  • 

uj  00  V)  . 

Sz  Jli 

</>  3 = 


ujDoa 
i—  _»  _i  o 
co  < <£  o 

in 

Ch  N ID  lO 
ID  rS  iD  ID 

r-4  CO  LO  LO 

ao<  o 


r—  __ 

« c 

v£> 

S—  O t/1 
O T-  0) 
O M-  +->  r- 

u, 

UJ 

CC  — J 

■X  JO  0 

0)  J-  3C 

E +■> 

-J-  01  -0 
(—  C QJ 

<U  x: 
uia  u 
> c 
>>  3 

(_>  u)  +->  Q. 

LlJ  3 •<-  01 

l/>  -r-  O S- 

X»  O a. 

JO  1— 

•>  CC  0)  JC 

UJ  > Ol 

JE  -*  J-  3 

*— ■ oa)0 

h-  o a.  1. 

l/l  X )- 


o co  <x> 


(ww)  sniavy  mohs 


30 


PRESSURE  AT  SHOCK  FRONT  (U,)  KILOBARS 


AFFDL-TR-75-102 


DIAMETER, 

PROJECTILE 

SHAPE 

VELOCITY, 

ENER 

MM 

MATERIAL 

KM/SEC 

J 

45.56 

STEEL 

CYLINDER 

4.27 

82' 

0 5.56 

COPPER 

CYLINDER 

3.96 

76; 

> 3.175 

STEEL 

SPHERE 

4.27 

11 

0 5.56 

ALUMINUM 

SPHERE 

2.32 

61 

a 1.59 

STEEL 

SPHERE 

5.35 

2: 

o 5.56 

ALUMINUM 

SPHERE 

1.90 

4' 

2 4 6 8 10  12 

SHOCK  RADIUS  CM 


Figure  10.  Shock  Front  Pressures  in  Water  as  a 
Function  of  Shock  Radius 
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Figure  11.  Shock  Front  Pressures  in  Water  (20°C) 
as  a Function  of  Shock  Front  Velocity 
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deduced  from  velocity  measurements.  The  energy  loss  in  the  fluid  due  to 
drag  can  then  be  calculated  reasonably  well  with  the  analytical  equation 
(6:19) 


(45) 


where  E0  is  the  energy  of  the  projectile  just  after  impact,  Cq  is  the 
coefficient  of  drag  of  the  fluid  on  the  projectile,/Df  and  pp  are  the 
densities  of  the  fluid  and  the  projectile,  respectively;  v0  is  the 
projectile  velocity  just  after  impact,  and  rp  is  the  projectile  radius. 
Equation  45  will  be  derived  later  in  a slightly  different  form  (Section 
III).  The  energy  loss  for  the  light  particles  studied  here  was  well 
over  90%  within  the  first  20  usee.  As  an  example,  an  aluminum  sphere 
with  3.175  mm  diameter  and  having  an  impact  velocity  of  5.27  mm/usec 
was  essentially  stopped  within  20  mm.  Thus,  it  is  natural  to  want  to 
approximate  the  energy  release  that  formed  the  shock  wave  as  a point 
source  for  these  very  light  projectiles.  Bullets,  even  fragments,  would 
not  normally  be  expected  to  follow  such  a large  energy  loss  percentage 
so  quickly. 

2.  COMPRESSION  OF  THE  FLUID  NEAR  THE  IMPACT  POINT 

Williams  (9)  has  presented  a model  of  the  fluid  compression  behind 
the  shock  front  based  on  the  fluid  volume  displaced  by  the  projectile 
during  the  first  couple  of  millimeters.  The  following  assumptions  are 
required:  (1)  the  projectile  penetrates  the  wall  smoothly  and 
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continuously,  (2)  the  projectile  and  shock  velocities  remain  constant 
during  penetration,  and  (3)  the  tank  wall,  volume  of  the  wall  removed, 
and  the  tank  wall  displacement  have  no  effect  on  the  fluid  behind  the 
shock  surface  (even  though  the  tank  wall  displacement  velocity  may  be 
significant).  With  these  assumptions,  then,  and  reference  to  Figure  12, 
the  total  fuel  volume  displaced  by  the  projectile  is  given  by 


vp  * ir /V(£) 

o 


(46) 


where  r(£)  is  the  distance  the  projectile  protrudes  into  the  fluid, 
R(£)  is  the  radius  of  the  projectile,  and£  is  a dummy  variable  which 
is  equal  to  0 at  the  tip  of  the  projectile  and  increases  to  the  left. 
The  volume  enclosed  by  the  hemispherical  shock  surface  is 


Vo  s 


2t r 


u;  r 


(47) 


The  ratio  of  the  volume  within  the  shocked  region  at  any  instant  and  the 
vo^rne  originally  occupied  by  that  region  is  the  compression 


I - 


(48) 
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Figure  12.  Projectile  - Shock  Front  Diagram 


where 


±.  . 4.  f*  '*({)«€ 

X,  T J»  uj  ts 


is  the  relative  mean  compression  ratio  of  the  fuel  within  the  shock.  In 
the  approximation  that  the  nose  of  the  projectile  is  a cone  with  angle 
9 , then  r(£)s£ton0  and  the  integral  is  easily  evaluated  to  give  the 
relative  mean  compression  ratio  for  this  case 


ton2  9 
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Thus  the  initial  mean  compression  ratio,  up  to  penetration  of  the 
shoulder  of  the  projectile,  is  essentially  time-independent.  This  model 
is  particularly  useful  for  the  calculation  of  very  early  shock  pressures. 

3.  DENSITY,  PARTICLE  VELOCITY,  PRESSURE  PROFILES,  AND  THE  ENERGY 
EQUATION 

If  the  projectile  were  to  release  all  of  its  energy  at  a single 
point  instead  of  continuously  along  its  trajectory,  one  could  reference 
the  published  works  of  the  authors  on  blast  waves.  This  is  basically 
the  limiting  case  that  Yurkovich  describes,  drawing  upon  the  work  of 
Bach  and  Lee  (5).  This  assumption  is  clearly  valid  for  light  particles, 
which  may  be  essentially  stopped  within  a few  centimeters.  However, 
this  model  can  be  criticized  since  a bullet  will  not  normally  lose  a 
large  percentage  of  its  energy  in  such  a short  distance.  This  is  es- 
pecially true  for  normal  impact,  which  is  the  case  generally  found  in 

testing  situations.  But,  it  may  not  be  a valid  criticism  for  tumbled 

* 

entry,  because  large  percentages  of  the  projectile  energy  could  be  lost 
during  impact.  And,  tumbled  entry  is  the  condition  which  can  be  most 
frequently  expected  in  a combat  situation,  where  the  bullet  must 
first  penetrate  the  aircraft  skin  at  an  oblique  angle,  and  with  a 
relative  tangential  velocity. 
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In  the  Yurkovich  model,  the  projectile  becomes  a point  source  of 
energy  immediately  inside  the  front  wall.  The  hemispherical  shock 
surface  that  Is  formed  is  taken  to  be  half  the  shock  sphere  formed  by  a 
blast  wave  in  an  infinite  fluid.  The  functional  form  of  the  density 
behind  the  shock  surface  must  be  known  or  assumed.  Then,  it  is  possible 
to  derive  profile  equations  for  the  fluid  density,  the  fluid  particle 
velocity,  and  the  fluid  pressure  in  the  shocked  region.  Finally,  the 
energy  in  this  region  can  be  written  in  terms  of  an  integral  which  must 
be  evaluated  numerically. 

a.  Density  Profile.  Consider  the  projectile  tip  to  have  just 
penetrated  the  front  wall  without  disturbing  the  fluid.  At  this  in- 
stant, the  projectile  deposits  a significant  amount  of  its  energy, 
resulting  in  a hemispherical  shock  surface  which  expands  radially  from 
the  impact  point.  This  hemisphere  will  have  the  same  characteristics  as 
half  a sphere  formed  by  a point  release  of  energy  into  an  infinite 
fluid;  i.e.,  all  wall  effects  are  ignored. 


An  expanding  shock  sphere  in  an  infinite  fluid  would  enclose  a mass 


m = 4ir 


P dr 


(51) 


where  p-p{ r,t)  and  r is  the  radius  of  the  sphere.  The  mass  which  had 
originally  occupied  this  volume  was 


m, 


3"  " 


(52) 
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But.  these  masses  must  be  equal,  thus 


r«  ,s 

/ '2/»dr  = ~ 


riP  o 
3 


or,  in  nondimensional  units  (Equations  30  and  33), 


/ / s (T,  Ms)df  * -y 


where  the  integration  can  be  performed  once  the  functional  form  of 
£ is  known.  Yurkovich  used  the  same  form  assumed  by  Porzel  (11) 


P = ?,(M*)7  Q 


Integration  of  Equation  55  together  with  the  result  given  in  Equation  54 
determine  q(Ms)  to  be  , 


q (M,)  » 3 </>,-!> 


The  density  profile  behind  the  shock  front  can  now  be  written  explicitly: 


3U-1) 

( r,  M,)  * py  ( Mg)  r 1 
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b.  Parti-le  Velocity  Profile.  Since  p is  known,  it  is  possible 
to  solve  the  Conservation  of  Mass  equation  in  nondimensional  coordinates 
foru  , the  nondimensional  particle  velocity.  Equation  38  is  rewritten 

+ pZf  * 2¥  8 (58) 

Substituting 


P 


S 


Py 


_q 

r 


dp_ 

dr 

dp 

dM, 


Jq  - 1) 
P , d r 


( *P. 

\aMfi 


(59) 


and  rearranging  terms, 


£1  + d+2  - 
dr  f 


q - /SM 


( 

8 V dM, 


Px  + 


in  r 


ag  \ 

dMs  / 


(60) 


Once  M#  is  fixed, q{Ms)  will  be  constant,  so  that  for  eachM#,  Equation 
60  can  be  rewritten  as  a first-order  linear  differential  equation  in  the 
standard  form 


^ + P (x)  y * Q (x) 


(61) 
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Using  the  integrating  factor  the  Integration  of  Equation  60  gives 


u _ 
7 


- /3m,  -i- 
9 dut 

1 

'f 

l 

q+3 
T 1 n7 

1 1 f C 

l q+3  ( 

q + 3)2-*  r«+3 

where  the  integration  constant  C is  found  to  be  0 for  the  boundary 
conditions  usu,  , and  ?=l  . Recalling  that  q+3»3/>,  . the  dertv-  1 

ative  can  be  written  as  ~ * 3 , and  the  particle  velocity  profile 

vm$  uM| 

has  the  final  form 


where 


u = U,  7 ( I - 0 I n7 ) 


“i  3m, 


Equation  63  can  be  simplified  for  strong  shocks.  Since/),  is  a 
function  of  M8  alone,  and,  further,  is  nearly  constant  for  shocks  which 
have  a constant  density  ratio  across  the  shock  wave,  the  velocity  ratio 
for  very  strong  shocks  is 

U a u.  7 (65) 
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c.  Fluid  Pressure  Profile.  For  both  density  and  velocity  pro- 
files, Yurkovich  has  followed  Bach  and  Lee  (5).  However,  in  obtaining 
the  pressure  profile,  he  makes  the  strong  shock  simplification  before 
the  derivation,  while  Bach  and  Lee  make  the  assumption  after  the  de- 
rivation. Since  the  final  equations  are  different,  and  since  neither 
has  been  verified  for  bullets,  both  will  be  presented  here  for  com- 
pleteness. 


The  conservation  of  momentum  is  restated  in  nondimensional  coor- 


dinates (Equation  43): 


To  find  the  normalized  pressure  functionP,  this  equation  is  rearranged 


for  integration  with  respect  tor  : 


P‘-fp[(*-7)  -|2 - + + " 


Noting  that/®,u  ,q  , and®  are  functions  of  Mt  alone,  the  following 

~ — q 

substitutions  are  made  P ~ P\f 

u = u,  7 (l  - 0 I n 7 ) 
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Rearranging  and  identifying  similar  terms,  gives 


P * - p J*[  P2  7**  + l + P37q+*  Inr  + P47fl  I n2r]  dr  (69) 


where  the  following  functions  have  been  defined: 


u,  (11,-1)  (1-6)  + £(u,  + M, 

rr,®  [i  ♦ u,(®-2  )-/3]-/3ms-^- u,® 

S 


IF,2®2 


and  P,  will  be  defined  later.  Equation  69  is  then  integrated  to  give 


— - / p2  JO +2)  Ps  (q  + 2)r  . _ l 

P = + (T^j2  l(l+z)lnr-lj 


♦ (7^',+2’([(-2)'—r+')}+ 


and  the  constant  of  integration  K,  is  determined  from  the  boundary 
condition  that  P=  P,  at  ni  : 


P<  ‘Pi 


f \ h , gp4  \ 
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Finally,  the  pressure  profile  can  be  written  as 


P = P,  + 


P,  - 


q+2  P2 


1 P, 


(<+2) 


2 ** 


{ 


■ + * 21 


[(  q +2  )ln 7 -l]} 


+ P4  { 2 - 7<q  + 2)  ( [ (q  + *)  In7  - I ]*+l)} 


(73) 


For  strong  shock  waves, p.s  — is  very  large  and  thus  u,  = l--»- 

P 0 P\ 

is  just  slightly  less  than  1.  The  third  and  fourth  terms  of  Equation 

73,  which  are  modified  by  terms  ,~~r:n  = -r.n  with  n=  2,  and  3, 

respectively,  can  be  eliminated,  and  P2  simplifies.  Thus,  for  strong 

shocks,  Bach  and  Lee  write 


P « P,  + ^ + 


(74) 


As  mentioned  earlier,  Yurkovich  makes  the  simplifying  assumption  of 
a strong  shock  wave  before  deriving  the  pressure  profile,  i.e.. 


p 

s 

P, rq 

c 

a 

V 

du 

dr 

a 

du 

at 

_ du, 

dMs 

at 

r dMt 

(75) 
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Putting  these  substitutions  into  Equation  67,  P can  be  written  as 


where  p£  is  just  Pg  as  in  Equation  70  when®  is  near  0: 


p*a  u,  ( u,  +)8- 1 ) +/3m* 


To  compare  Yurkovich's  results  with  those  of  Bach  and  Lee  (Equation 
74),  Equation  76  is  written  out  as 


♦*-•>♦  (78) 

' 1 q + 2) 

The  last  term  in  Equation  78  is  an  extra  term  for  Yurkovich's  results; 
no  study  of  its  significance  was  made  for  this  report. 


d.  Energy  Equation.  A final  relationship  that  is  available  for 
describing  the  properties  within  the  shocked  region  is  the  energy  equa- 
tion. Since  the  total  energy  deposited  into  the  fluid  is  conserved,  the 
energy  deposited  in  the  hemisphere  is  the  sum  of  the  change  in  the 
internal  energy  and  the  change  in  the  kinetic  energy.  For  a sphere, 

E»s  2*/  [(•-•<>)/>  ♦ U "2 U°  /»]  r*dr  (79) 
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or,  for  a hemisphere  in  nondimensional  coordinates  with  u0=  0 , 


Eh.,* 


*7>oCo  M?  '•*/[(*•  + "1-  ] "r 


(80) 


where  the  only  quantity  not  known  within  the  integral  is  the  specific 
internal  energy 5.  This  work  is  tremendously  simplified  for  the  strong 
shock  assumption;  I0  is  small  and  can  be  neglected.  Further,  for  a 
strong  shock,  the  fluid  variables  will  remain  relatively  constant  for 

some  distance  behind  the  shock  surface,  thus  permitting®  to  be  replaced 

-2 

bye,.  Combining  appropriate  relationships.  Equation  41  (i,**^*)  , 

Equation  55  [P-P,fQ)  , and  Equation  65  (u*u,r ) , to  rewrite  Equation 

80,  gives 


Eh.s  = f - />0C © Mf  r s “*  P\  £ ( I *' r* > rfa+2)df 


(81) 


Integrating  over  the  hemisphere,  and  replacing  q with  3(/>,-l)  , 
finally  yields 


Eh.*.s  '3'C*  r*  p0  u* 


(2fl±L\ 

' 3A>+2  l 


(82) 


e.  Discussion.  In  summary,  profiles  have  been  developed  for  the 
density,  particle  velocity,  pressure,  and  energy  of  the  fluid,  the 


45 


•I.’.-uAr...,;*?' 


AFFDL-TR-75-1 02 


pressure  and  energy  relationships  being  the  ones  which  will  be  applied 
most  often.  But,  p,  ,/>,  , u,  , and/3,  have  not  yet  been  evaluated. 

Yurkovich  has  used  the  Tait  equation  of  state.  Equation  42;  so 
that  the  value  of  P,  just  behind  the  shock  surface  can  be  found  from 

[;h  <•» 


which  from  Equation  40  must  also  be  the  expression  for  u,  . Now  Equation 
39  cun  be  used  to  get  a relationship  for p,  : 


Pi 


I 


I 


n Ms 


(84) 


A direct  solution  of  Equation  84  for  pt  in  terms  of  n and  M«  ns  not  pos- 
sible, however,  a table  of  values  of  /»,  vs  Ma  for  a given  value  of  n 
can  be  computed  from  a rearrangement  of  Equation  84.  The  relation 


fi-i 


(85) 


was  used  to  plot  Figure  13.  Then,  from  Equation  39,  P,  and  uf  are 
easily  found  also  as  functions  of  n and  Ms  from 

P,=  U,=  l--Jr  (86) 

Pi 

Equation  86  is  plotted  in  Figure  14. 
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With  P,  known.  Equation  33  can  be  used  to  find  the  pressures  just 
behind  the  shock  front  from 


pt  c pt  Po»s8 


This  information  is  plotted  in  Figure  15. 


To  solve  for  the  pressure  profile,?  , it  is  necessary  to  find 


2 n Ms  </>,  - IT 

- h ♦ i ~~~  ^T5 
np(  py  +1 


a©  , jSm*  dV,  ^ f i , d$  t aP, 

Im7  77  W + ® (89) 

The  term#  is  taken  as  constant  by  Yurkovich,  and  with  this  assumption, 

the  term  involving  is  zero  in  Equation  89. 

<?nn# 


Finally,  a value  for  # can  be  estimated  from  Equation  38.  Based 
on  Figure  8,  the  shock  radius  has  already  been  assumed  proportional  to 
time  to  some  power  a , which  has  been  approximated  by  0.8  for  light 
hypervelocity  particles.  Thus,  by  assuming  r*»hta  , Equation  38  gives 


£ * >-± 


and,  for  meteorite  work  at  least,  # = -0.25. 
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Figure  15.  Pressure,  P-j  at  Shock  Front  as  a Function  of 
the  Shock  Mach  Number,  Ms 

The  pressure  function  P can  now  be  evaluated  from  Equation  76  and 
plotted  as  in  Figure  16.  These  curves  show  the  pressure  profiles  based  on 
the  strong  shock  assumption  as  implemented  by  Yurkovich.  Two  other 
pressure  profiles  are  available  from  the  work  of  Bach  and  Lee:  a com- 

plete solution  (Equation  73)  and  another  solution  based  on  a strong 
shock  assumption  (Equation  74).  These  three  solutions  are  plotted  in 
Figure  17  for  n = 7,  (3  * -0.25,  and  Ms  = 1.5.  The  curves  for  the  other 
mach  numbers  between  1 and  2 have  also  been  plotted,  though  they  are  not 
pictured  here.  They  all  show  the  same  forms  as  in  Figure  17,  and  with 
lower  pressures  for  lower  shock  mach  numbers.  Only  the  Bach  and  Lee 
strong  shock  curve  does  not  become  negative  for  weaker  shocks,  indi- 
cating that  an  adjustment  is  required  before  these  relations  are  ap- 
plied, especially  to  cases  of  weaker  shocks,  as  with  projectiles. 
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Figure  16.  The  Nondlmenslonal  Pressure,  F (Ms,  r),  as  a Function 
of  Nondimensional  Distance  for  Several  Values  of 
Ms  and  for  n = 7 

One  possible  critical  assumption  was  in  the  evaluation  of  the  term 
0 . The  term  0 was  treated  as  a constant;  and  the  measurements  of  the 
shock  front  radius  (Figure  8)  do  indicate  that  a (hence  0 ) is  a con- 
stant. But  these  curves  were  obtained  for  light  hypervelocity  projec- 
tiles for  a very  short  time  interval  and  there  is  no  reason  to  assume 
that  these  curves  will  be  so  nearly  straight  for  either  the  lower  shock 
front  velocities  of  bullets  or  for  the  longer  time  intervals  that  may  be 
Important  for  fragments.  The  solution  suggested  by  Yurkovich  did  not 
present  the  method  for  dealing  with  a varying  0 , and  so  is  not  dis- 
cussed in  this  report.  But  j8  Is  a function  of  Mt  and  the  technique  for 
finding  /3(MS)  is  available  in  Bach  and  Lee  (Reference  5).  It  must  be 
recommended  therefore,  that  before  applying  this  analysis  to  the  shock 
portion  of  the  hydrodynamic  ram  effect,  either  a should  be  measured  to 
determine  its  value,  or  preferably  that  this  analysis  should  be  extended 
to  Include  the  method  used  by  Bach  and  Lee  to  find  the  functional  form 
of  0 . 50 
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| Another  result  of  Yurkovich's  analysis  of  the  shock  phase  is  the 

energy  density  equation  which  Is  gained  by  dividing  Eht  in  Equation 
81  by  the  volume  of  the  hemisphere  enclosed  by  the  shock. 


Equation  90  is  plotted  In  Figure  18  for  the  condition  n 3 7 (water)  and 
gives  the  energy  density  at  the  shock  front.  Since  the  derivation  of 
Equation  87  was  based  on  the  strong  shock  assumption,  the  dependence  of 
r is  lost  from  the  relationship,  and  hence  there  is  no  energy  profile. 
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SECTION  IV 

CAVITY  DESCRIPTION  AND  FLUID  DRAG  ANALYSIS 

No  aspect  of  the  hydrodynamic  ram  phenomenon  has  received  more 
attention  than  the  cavity  phase.  The  motion  of  the  fluid  has  been 
monitored  inside  the  tank  du. ing  thousands  of  tests  by  high  speed 
cameras  (seldom  over  7000  fps)  peering  through  plexiglass  windows, 
and  by  pressure  transducers  recording  high  pressures  and  sometimes 
shocks  at  half  a dozen  or  more  points  throughout  the  tank.  Strain 
gages,  witness  sheets,  and  even  cavity  gages  have  been  used  to  moni- 
tor the  violent  fluid  and  tank  motions  for  up  to  a tenth  of  a second. 

As  the  projectile  penetrates  the  tank,  the  bullet  displaces  the 
fluid  by  imparting  velocities  to  the  fluid  particles  in  the 
direction  of  the  normal  to  the  bullet  surface  at  every  point.  The 
observed  net  result  is  that  the  fluid  is  forced  to  move  radially 
away  from  the  bullet's  path,  leaving  a cavity  in  the  fluid  and  creat- 
ing high  pressure  in  the  vicinity  of  the  bullet  of  35  to  140  bar 
(400-2000  psi)  at  about  10-20  cm  from  the  bullet. 

Occasionally,  the  projectile  will  remain  stable  throughout  its 
trajectory  and  the  cavity  that  is  produced  is  just  a cylinder  about 
the  path.  When  this  happens,  there  is  little  visible  motion  of  the 
fluid  and  the  projectile  loses  a minimum  amount  of  velocity  before 
it  exits  the  tank. 
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Normally,  however,  the  bullet  will  tumble  within  20  to  30  cm.  The 
larger  the  drag  the  fluid  exerts  on  the  bullet,  the  more  projectile 
kinetic  energy  is  transferred  to  fluid  kinetic  energy.  Thus  a 
tumbling  0.50  caliber  bullet  can  be  stopped  by  a 3-foot  tank,  while 
a stable  projectile  may  not.  As  the  tumbling  projectile  gives  up  vary- 
ing amounts  of  its  energy,  the  diameter  of  the  cavity  it  creates 
changes  proportionally;  a smoothly  tumbling  projectile  will  create 
a series  of  connected  cavities  with  regularly  varying  cross  sections, 
with  the  location  of  the  smallest  cross  sections  corresponding  to 
the  points  in  the  trajectory  where  the  bullet  is  traveling  nose 
first,  and  the  largest  cross  sections  corresponding  to  those  parts 
of  the  trajectory  where  the  bullet  is  traveling  with  its  length  fac- 
ing forward.  The  cavity  grows  into  a sphere  as  the  smaller  cross 
sections  begin  to  collapse  while  the  wider  cross  sections,  whose 
fluid  particles  have  more  kinetic  energy,  are  still  expanding.  At 
its  maximum  radius,  the  cavity  is  most  nearly  spherical  and  remains 
spherical  as  it  collapses  to  some  minimum  radius.  If  the  cavity  is 
fairly  large  and  nearly  spherical,  it  will  rebound  to  a second  maxi- 
mum radius  which  is  about  2/3  the  diameter  of  the  first  cavity. 

This  cycle  can  be  repeated  two,  three,  or  more  times  before  the 
cavity  dies.  Typically,  a 0.50  caliber  bullet  can  be  expected  to 
have  about  a 50  cm  first  maximum  radius  and  about  a 5 cm  first 
minimum  radius.  The  first  minimum  can  be  expected  at  about  30  msec 
after  impact,  and  a second  minimum  about  30  msec  after  the  first. 
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In  this  section  of  the  report,  methods  are  described  for  estimating 
the  pressures  in  the  fluid  from  a knowledge  of  either  the  cavity  radius 
or  the  cavity  period.  Equations  which  relate  the  fluid  energy  to  the 
projectile  energy  are  developed,  and  equations  which  describe  the  pres- 
sures in  the  fuel  tank  as  functions  of  time  and  position  due  to  pro- 
jectile drag  are  derived.  Following  these  are  some  brief  remarks  on  the 
pressures  associated  with  cavity  collapse  and  then  a concluding  dis- 
cussion. 

1.  PRESSURE  ESTIMATION  TECHNIQUE 

Cardea  and  Torvik  (Reference  12)  present  a method  of  estimating  the 
pressure  in  the  tank  as  a function  of  the  distance  from  the  center  of 
the  cavity  if  either  the  maximum  diameter  of  the  cavity  om  or  the  cavity 
period  is  known.  Much  of  their  work  is  based  on  Cole  (3:364),  who  gives 
the  maximum  recorded  pressure  gage  reading  in  terms  of  the  weight  of  TNT 
exploded  by 


Pgm 


2590  wi/s 
r 


(92) 


where  W is  the  weight  of  the  TNT  in  pounds, r is  the  distance  from  the 
detonation  in  feet,  and  Pgm  = pnTpo  in  Psi-  While  Cardea  and  Torvik 
note  that  the  energy  in  one  pound  of  TNT  is  \ .514  x 106  ft-lbs.  Cole 
observes  that  about  half  this  energy  goes  into  the  formation  of  the 
cavity.  This  leads  to 

1931 

where  Ec  is  the  energy  required  to  form  the  cavity. 
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There  are  two  methods  to  estimate  Ec.  By  assuming  that  the  cavity 
behind  the  projectile  is  sufficiently  like  that  produced  by  an  under- 
water explosion  of  TNT,  Cole's  results  can  be  used.  Thus,  for  a cavity 
in  an  incompressible  fluid  sufficiently  distant  from  the  tank  wall  to 
ignore  reflection  the  energy  required  to  form  the  cavity  is  a function 
of  the  maximum  cavity  radius om,  and  is  given  by 

Ec  ■ -if-  Po  °n?  (94) 


where  P0  is  the  ambient  pressure.  For  use  in  Equation  92,  Equation  93 
is  rewritten  as 


Ec,/8  = 1.612  om  P0I/8 


(95) 


A second  method  for  estimating  the  cavity  energy  is  given  also  by 
Cole  (3:276)  as  the  Willis  formula 


p #/• 

no 

U4poUt 


T 


(96) 
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where  T Is  the  time  for  one  cavity  oscillation,  and/so  Is  the  ambient 
density.  This  second  method  Is  preferable  because  the  cavity  period  Is 
measured  more  easily  than  the  maximum  cavity  radius. 


2.  PROJECTILE/FLUID  ENERGY  BALANCE 

As  the  projectile  moves  through  the  tank,  its  motion  is  opposed  by 
fluid  particles;  the  total  force  opposing  the  projectile  is  the  fluid 
drag.  It  is  assumed  that  all  of  the  energy  that  is  given  up  by  the 
bullet  is  either  absorbed  by  the  fluid  particles  as  fluid  kinetic  energy 
or  lost  as  work  done  on  the  fluid  in  displacing  the  particles  and 
creating  the  cavity.  Thus,  the  equation  to  solve  is 


dEp  = dEf  + dW 


where  the  work  done  against  external  pressure  Pe  to  change  the  volume  i: 


dW  s P#  d V 


Considering  the  volume  expansion  along  each  increment  of  length  dx  to  be 
always  perpendicular  to  x , 


dW  * (P0-  Pc)  7 r o*  dv 
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Both  Yurkovich  (8)  and  Lundstrom  (13)  have  completed  the  solution 
of  Equation  96  In  detail.  There  are  fundamental  differences  to  their 
approaches,  but  their  basic  schemes  are  parallel,  and  they  both  derive 
similar  incidental  information  in  the  process.  To  illustrate  the  simi- 
larities of  their  approaches  and  to  facilitate  a comparison  of  the 
results,  their  solutions  for  each  of  the  major  steps  of  the  energy 
balance  equation  will  be  given  sequentially.  Thus,  both  solutions  of 
dEp  are  given,  and  then  both  solutions  of  dEf  are  given.  Preceding  the 
calculation  of  dEf  , it  is  necessary  to  show  the  derivation  of  the  source 
potential  <f>. 

a.  Calculations  ofdEp.  The  kinetic  energy  lost  by  the  projec- 
tile is  given  by 


dEp  a - mv  dv 


(100) 


i 


■4 


wherem  and  v are  the  projectile  mass  and  velocity,  respectively.  Then 
the  energy  loss  along  the  trajectory  is 


dEp 

dx 


- mv 


dv 

dx 


S 


(101) 


But  the  last  term  is  the  force  opposing  the  projectile,  i.e.,  the  fluid 
drag 


Sr  ■ 0 ■ 


Ap  (x)  Cq  (x) 


(102) 
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where  p is  the  fluid  density,  CD(x)  Is  the  coefficient  of  drag,  and 
Ap(x)  is  the  area  that  the  projectile  presents  to  the  fluid  as  it  moves 
along  its  trajectory.  Both  the  presented  area  and  the  coefficient  of 
drag  change  as  the  projectile  moves  through  the  fluid;  however,  the 
problem  Is  normally  simplified  by  fixing  the  presented  area  of  the 
projectile  and  considering  only  the  coefficient  of  drag  as  a variable. 

The  presented  area,Ap(x)  , of  the  projectile  is  taken  to  be  the  area  of  its 
widest  cross  section,  Ap . 

Lundstrom  combines  Equation  101  and  Equation  102  to  get 

V - -£dx  (103) 


where  the  velocity  decay  coefficient  Is  defined  as 

s • j?ArSp 


(104) 


Equation  103  is  integrated  to  obtain v , 


v0  • 


-/£  d> 


(105) 
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Finally,  by  combining  Equations  101  and  104,  Lundstrom's  final  equation 
for  the  energy  deposition  along  the  trajectory,  can  be  written  as 


mv*  0 


(106) 


Equation  105  can  be  used  to  write  the  final  equation  (Equation  107)  In 
terms  of  the  constant  v0  . which  is  measurable  at  Impact,  and  the  func- 
tion 0 which  must  be  determined  by  observation: 


dEp  -«/£dx 

-dT  * mvo 


(107) 


Yurkovich  combined  Equations  101  and  102  differently  to  get 


^ - -/3dt 
v* 


(108) 


This  equation  is  not  integrated  to  findv  since  the  functional  form  of 
0 is  not  known.  Yurkovich  approximates  0 with  0\  , the  average  value 
of  0 in  the  ith  segment  of  the  projectile  path.  Then,  Equation  108  is 
between  the  end  points  of  each  segment  to  give 


i + £„.vi  At 


009) 
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where  the  subscripts  i and  l+l  refer  to  values  at  the  beginning  and  end 
of  the  1th  segment,  and  At  * t|+i -t  \ . 0j+|  Is  now  an  unknown  constant. 

Yurkovlch  uses  Equation  109  to  find  two  more  relationships. 

Equation  109  is  Integrated  to  find  the  distance  the  projectile  moves 
during  each  time  Interval, At  , 


Ax  ■ In  [ I + /3I+|  v,  At  ] 


0 


(HO) 


Ui 


Also,  Equation  109  can  be  combined  with  Ep*  to  get  the  energy  of 
the  projectile  in  each  segment 


Ep, 


Oil) 


(For  constant  Cg,  this  would  reduce  to  Equation  45,  which  was  given 
without  derivation.)  Finally,  the  total  energy  loss  due  to  drag  is  the 
sum  of  all  the  segments  of  energy 


Eo  - E 


n*i 

pf  ‘ £ E| 


012) 


Equation  112  can  be  differentiated  numerically  to  findY  , the  energy 
transferred  per  unit  length  of  the  projectile  path.  ThendEp  is  approx- 
imated by 


AEd  * YAx 


(113) 
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b.  Velocity  Potential  Due  to  a Line  of  Sources.  The  final  goal 
of  the  fluid  drag  section  is  the  solution  of  Bernoulli's  equation. 
Equation  28,  for  the  pressure  distribution  throughout  the  fuel  tank,  and 
this  equation  requires  the  differentiation  of  the  velocity  potential, 

<f>.  Thus  both  authors  have  described  the  fluid  energy  differential, 
dE#,  in  terms  of^>.  Yurkovich  has  drawn  upon  the  known  velocity  po- 
tential of  a spherical  bubble,  while  Lundstrom  has  derived  an  appro- 
priate velocity  potential  for  the  problem  based  on  the  concept  of  source 
strengths.  Lundstrom's  velocity  potential  is  derived,  in  the  following 
paragraphs . 

.#■ 

A 

The  potential,^,  can  be  constructed  from  a line  of  sources  defined 
as  points  in  the  fluid  from  which  fluid  flow  can  be  thought  to  emanate. 
These  sources  are  imagined  to  exist  at  close  intervals  along  the  pro- 
jectile trajectory  as  in  Figure  18.  While  the  concept  is  purely  fic- 
titious, it  can  be  used  effectively  to  describe  the  flow  of  fluid  away 

\ 

from  an  explosion  or  from  the  trajectory  of  a projectile.  The  source  is 
assigned  a strength.m  , which  is  related  to  the  rate  of  fluid  flow 
(msAv/At)  • The  fluid  flow  associated  with  a point  source  is  de- 
termined by  that  volume  of  fluid  which  would  pass  through  the  surface  of 
an  imaginary  sphere  enclosing  the  source  (Figure  19): 

m At  « [<r  + Ar)*~r»]  (114) 


which,  after  simplifying  and  taking  the  limit  leaves 

m ■ 4 ir  r*  f 
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Since  u«1  the  flow  is  radial,  the  velocity  potential  gradient  must  be 
just  r , so  that 


f!£  . 

dr  4irr* 


(116) 


One  integration  shows  that  the  potential  function  for  the  point  source 
of  fluid  is  given  as 


* 


m 

47rr 


(117) 


The  fluid  velocity  away  from  the  trajectory  of  a penetrating  bullet 
can  be  thought  of  as  due  to  many  sources  along  the  bullet  path,  each 
having  the  above  potential.  Thus,  the  potential  to  solve  in  Bernoulli's 
equation  is  > for  arbitrarily  close  sources,  it  is 

♦ ■ /Xj2£Ld<  (m) 

0 


where  £ is  a dummy  variable  of  integration  along  the*  axis  andX 
is  the  interval  of  projectile  travel. 

A final  simplification  is  made.  Nearly  all  of  the  fluid  velocity 
is  observed  to  be  radially  outward  from  the  trajectory;  i.e.,  all  of  the 
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motion  Is  confined  to  the  slice  d£  for  each  Increment  of  the  bullet 
trajectory.  This  allows  the  source  strength  to  be  represented  by  a 
function  m's—f , which  represents  the  flow  rate  In  each  cross  section 
of  the  trajectory: 


m'  At  « tt  [<R  ARl*  -R*] 
In  the  limit,  m'  is  simplified  to 

m1  « 2irRR 


(119) 


(120) 


Equation  118  is  finally  reduced  to 

<f>  * d£ 

o 

where  a source  strength  function  $ was  defined  by 

C « RR 


(121) 


(122) 


c.  Calculations  of  dEf • Yurkovich  approximates  the  cavity  behind 
the  passing  projectile  as  a set  of  expanding  spheres,  with  centers  fixed 
on  the  projectile  path  (assume  the  x axis  for  simplicity).  Then  the 
velocity  potential  outside  this  expanding  gas  sphere  is  known,  and,  with 
reference  to  Figure  18,  the  energy  of  the  fluid  can  be  calculated. 
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The  kinetic  energy  of  the  fluid  due  to  the  expanding  bubble  follows 

from 

6f  ■ / £(?<*»' <f» 

(123) 

where,  for  radial  flow,  the  velocity  has  been  replaced  by7</>.  Lamb 
(14:46)  derives  the  relationship  that 

f <V^)*dv  * -f  ds  (124) 


and  also  gives  the  velocity  potential  (14:122)  as 


<*>  = 


a2  6 


(125) 


for  a bubble  with  radius  a . 
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It  is  possible  now  to  evaluate  the  fluid  energy  at  the  bubble 

surface.  At  r * a:  = - 6 and  dS * o*ln  dd9d (0  Thus 

dr  Ir  s o * 

jf  *-^-dS  = | | (ao)(-6)ofln0d0d/  (126) 

so  that 

Ef|r  = o * 027) 

Noting  that  the  volume  of  the  spherical  bubble  is  , Equation  127 

o 

can  be  written  as 


028) 


Yurkovich  assumes  that  this  relationship  will  hold  for  volumes  of  arbi- 
trary shape,  i.e.,  for  cylinders  of  length  A*,  thus 


AEf  = \PQ"  (oo)2  0 29) 


Now,  combining  Equations  97,  99,  113,  and  129,  Yurkovich  writes 


that 


aep  = ya*  * p0  (°  ®)2  A*  +(%  “ pc)  A* 
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where  is  included  here  for  comparison  with  Lundstrom's  expressions 
(Yurkovich  allowed  P.  to  be  zero). 

Lundstrom  follows  Birkhoff  (15:228)  in  expressing  the  fluid  kinetic 
energy  throughout  the  cross  section  dA  of  the  trajectory.  The  fluid  is 
considered  inf  inite  and  is  bounded  only  by  the  cavity  radius  a (£).  But 
since  all  the  flow  has  been  assumed  radial,  and  since  it  will  be  neces- 
sary to  assume  only  noncompressible  fluid,  an  upper  limit  of  the  integration 
017,  where  17  >|  , is  established.  Thus,  the  fluid  kinetic  energy  in  the 
slice  dA  is  given  by 


dA  I dx 


(131) 


Using dAs  2voda  and  evaluating  £ at  the  cavity  boundary  allows  ft2dA 
to  be  written  as  -*j—.  But  £ is  proportional  to  the  line  source 
strength,  which  in  turn  is  constant  for  each  value  ofxs£  . The  in- 
tegral in  Equation  130  can  now  be  evaluated  to  give 

dEf  s wp  £2  1017  dx  (132) 


Lundstrom  adds  that  for  17  in  the  20  to  30  range,  correct  cavity  shapes 
are  predicted;  thus,  1017  is  between  3.0  and  3.4  and  can  be  considered 
constant. 


rr.rcrxr::;.  "va*:!;..  / ir;:^ 
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Having  evaluated  the  kinetic  energy  in  each  slice,  it  is  now  possi- 
ble to  write  the  energy  transfer  rate.  Defining  y as  the  energy  trans- 

dEp 

ferred  per  unit  distance, Y * — ■-  , gives 

dx 

dEp  « V dx  (133) 


Then,  by  combining  this  with  Equations  132  and  98,  the  energy  balance 
equation  can  be  written  in  the  form 


dEp  a Ydx  = ir/>  £ In  77  dx  + ^ I*  - irtT  dx 


(134) 


Yurkovich's  simplified  approach  has  yielded  a difference  equation  which 
compares  very  well  to  Lundstrom's  differential  equation  except  for  the 
first  term  on  the  rig:.t.  For  Ini ? in  the  range  3.0  to  3.4,  the  term 
irp£  In  rj  in  Equation  134  is  more  than  twice  the  term  •^•/>(oo)2  in 
Equation  130. 

d.  Additional  Results.  The  energy  balance  equation  serves  two 
purposes.  First,  it  can  be  rearranged  and  differentiated  to  yield  ex- 
pressions for  *,he  cavity  expansion  rate  and  the  maximum  cavity  radius. 
Second,  in  Lundstrom's  derivation,  the  energy  balance  equation  must  also 
yield  a suitable  form  for  the  source  strength  function  which  can  be 
integrated  to  give  the  velocity  potential.  Only  the  differential  equa- 
tion, Equation  134,  will  be  used  in  deriving  the  following  expressions 
though  similar  expressions  could  be  obtained  from  the  difference  equa- 
tion, Equation  130. 
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Equation  134  can  be  solved  for  £ : 


C*  - 2 .*] 

P L (w  R -f>)  1 


'o  ’Cl 


(135) 


Then,  with  the  definitions 


m 


(136) 


and 


P - Pr 
0 c 

p ,n  v 


(137) 


the  source  strength  function  becomes 

£ * T ^ - oZ  (138) 


Later,  <*„,(£)  will  be  shown  to  be  the  maximum  cavity  radius  at  £ , and 
a can  be  thought  of  as  the  average  cavity  expansion  rate,  though 
neither  term  is  quite  constant  since  may  vary  slightly.  The  negative 
solution  was  ignored  for  Equation  138  because  the  sources  must  be  posi- 
tive for  cavity  expansion;  and,  cavity  collapse  cannot  be  described  by 
the  velocity  potential  given  in  Equation  121. 
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At  the  cavity  surface  (or  for  R between  o and 017  for  Incompressible 
flow),  £ can  be  written  as  £ = oa  from  Equation  122.  Making  this  sub- 
stitution in  Equation  133  and  rearranging  leads  to 


a 6 


(139) 


from  which  the  cavity  expansion  rate  can  be  written  as 


■ r - 


(140) 


Equation  139  can  also  be  used  to  derive  a relationship  for  the 
source  strength  function  £ which  does  not  depend  upon  the  cavity  radius 
a . Rearrangement  and  integration  of  Equation  139  from  Ip  , the  tire 
the  bullet  arrives  at  xp,  tot  , i.e.. 


(141) 


yields 


V- 


m 


- a' 


= am  - 


(•  -v) 


(142) 
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Thus,  the  source  strength  function  can  be  written  as 

C 8 [an  - 0 (t  * tp)J  <143* 


Since  the  maximum  cavity  radius  occurred  at  6*0,  It  Is  apparent 
from  Equation  139  that  the  maximum  radius  is  justam.  And,  from  Equa- 
tion 142  the  cavity  reaches  Its  first  maximum^#  « om)at  timet,*. 


tp  + 


(144) 


Rewriting  Equation  144  to  get 

' « ('n  ' rp)  (145) 


a is  seen  to  be  just  the  average  growth  rate  of  the  cavity  during  the 
first  expansion. 

3.  PRESSURE  DISTRIBUTION 

The  pressure  field  generated  in  a fuel  tank  by  a decelerating 
projectile  is  given  by  Bernoulli's  equation  (Equation  28)  for  the  ir- 
rotational  flow  of  an  incompressible  fluid  with  no  viscosity.  The 
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derived  expression  Is  abbreviated  by  neglecting  external  forces,  and  the 
constant  of  Integration  Is  determined  from  the  Initial  static  condition 

p 

to  be  F 1 1 ) —-2-.  The  equation  to  solve  Is 
P 

P-P0  = _ (7<ft)2  _ dj>  (146) 

p 2 di 

Both  Yurkovich  and  Lundstrom  solve  this  equation,  Lundstrom  uses  the 
velocity  potential  derived  above  while  Yurkovich  (3:289)  uses 

* * - -4*-  - 4-  *"•  047) 


where  9 is  the  angle  between  x and  r as  shown  in  Figure  22. 

In  the  solution  of  the  pressure  distribution,  using  Equation  147, 
the  term is  given  in  spherical  coordinates  by 

‘ -5T  + + if-  (,48) 

In  the  earlier  derivation  the  center  of  the  sphere  was  fixed,  but  for 
this  portion  of  the  problem  the  sphere  must  be  allowed  to  move  in  the 
fixed  tank.  Therefore  the  term— 2- must  be  written  as 


5 (J+-\  . „ (J£S 

di  \ di  )m  * \ 3x  /t 
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Figure  22.  The  Angle,  0,  for  a Sphere  with  Forward  Motion 


where  m 
forward 


refers  to  the  fixed  coordinate  system  of  the  tank  and  »#  is  the 
velocity  of  the  sphere.  Finally,  it  is  necessary  to  write 
in  spherical  coordinates: 


cos 


d(t>  _ sine 

dr  r 


J± 

d6 


(150) 


The  following  then  is  Yurkovich's  solution  of  Bernoulli's  equation  for 
an  expanding  sphere  moving  horizontally  in  the  fluid  and  at  some  dis- 
tance from  the  tank  walls: 


~ 8 + dV  (°2®)  + i (+)*  co»®  (°*s  + 5v8o) 

Po 

+ ("H*  vs  (cos2e  - i sin2e)  - (-f-)4  (151) 

“ (“r)9vs  o co»®  - (-?“)*  -o'  (co*2®+  fsin2e) 
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Yurkovlch  notes  that  this  equation  can  be  used  not  only  for  cavity 
growth,  but  also  for  cavity  collapse,  provided  great  care  is  used  In 
choosing  the  proper  parameters  for  the  velocity  potential. 

Equation  151  can  be  simplified  for  two  important  cases  for  which 
v#  is  small.  First,  for  distances  far  from  the  cavity,-^  and  higher 
order  terms  are  sufficiently  near  zero  that 

S 4“  "at  ( o*  Q 1 + l ■f"  )*  co*  8 (152) 

P o 


If  in  addition,  there  is  no  forward  velocity, * 0 , and 


P-Pq 

Po 


a T ^ ta*6  1 


(153) 


This  final  condition,  lack  of  forward  velocity,  is  characteristic  of 
cavity  collapse. 

Lundstrom's  solution  of  Bernoulli's  equation  is  based  on  the  ve- 
locity potential  derived  earlier  from  a line  of  sources,  rewritten  as 

T ,xp<r)  om<e>-5rt-4--tb<e>l 
<f>  ( x,R,t,)  ■ -f  / — -===== -d£  (154) 

T'°  y-i'+t.-fi* 


where  the  terms  are  defined  as  in  Figure  23.  The  form  for  the  source 
strength  function  £ is  taken  from  Equation  143;  om(£)and  r(£)are  func- 

75 


•vgt 


;•*  wr.  *■ 


I:  I 


/a 


a , 
/ 

/ 

't  y 


/|  ( xi(°) 

/ 

/ 

. y ' 

. / FLUID 

S 

CAVITY 


5-Xb(t) 


Figure  23.  Flow  Diagram  for  the  Drag  Phase 

tions  of  position  along  the  projectile  path;  andXb(v)is  the  position 
of  the  projectile  at  the  retarded  time  r « t - r/c.  The  partial  de- 
rivatives of  ^ can  be  calculated  by  using  Liebnitz's  rule  (16:274)  which 


states  that  if 


,b(t) 

F(  r,  t ) «/  f (£,♦> 


055) 


i «M,ir 


(156) 
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Using  Llebnitz's  rule  on  Equation  154  yields 


a*  Z rap  axp(r)  T/.Xp(r) 

7T  “”Tl^‘“ar"“0/ 


ie]  (157) 


where  QpaOm^-xp  and  RpsR*+[)l“^pw]*  . The  Integral  In 

Equation  157  can  be  evaluated  from  tables  as  tn [ , -^p»Rp]  where  Ro 
is  the  distance  from  the  point  (x,R)  to  the  impact  point.  It  only 

remains  to  find  a suitable  expression  for  the  term  4~ p.  Noting  that 

0* 

f^Pav  , the  bullet  velocity  evaluated  atXp,  gives 

OT  P 


dip 

~JT 


* V, 


-dr_ 

p at 


(158) 


But  r «♦  - r/c  leads  to 


I'- 


4*-  ■ i + -t 

at  c 


axP 

dt 


059) 


dr 


Finally,  taking  jrl  __  . . gives 

0 1 X =ip(T) 


II 


rl- 

R 

■'  * yi 

£> 


dXp(T} 

dT~ 


I - 
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Then  the  final  form  for 


Is 

d* 


A + 


In 


x - Rj 

Xp(r)  + Rp 


where  A is  defined  as 


(161) 


A * 


X - Xp(r)  Vp 
jjjjT  •“  - 1 


(162) 


The  derivatives 


are  analogous. 


The  results  are: 


d£ 

dt 


A 


x - Xp(T) 
2e  Rp 


,XP(r) 

+ if  C (0 


(163) 


and 


- a ■ — 8 — + _a_  /Ip(Tl 

dR  2c  Rp  2 JQ 


(164) 
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where  £ Is  defined  as  In  Equation  143.  The  integrals  cannot  be  eval- 
uated until  the  form  of  £ is  known,  and  £ cannot  be  predicted  unless 
some  function  is  chosen  for  the  tumble  of  the  projectile. 


Equations  161,  163,  and  164  can  be  used  to  rewrite  Bernoulli's 
equation  (Equation  146): 

4.  CAVITY  COLLAPSE  PRESSURES 

Collapse  of  the  cavity  is  the  final  phenomenon  associated  with 
fluid  drag.  A 50-cal  projectile  will  exit  a one  meter  tank  in  less  than 
2 msec;  the  formation  of  the  cavity  is  still  taking  place  at  this  time, 
and  collapse  does  not  occur  until  about  30  msec.  By  the  time  collapse 
begins,  the  cavity  normally  has  a nearly  spherical  shape  and  very  little 
forward  motion.  Thus  the  fluid  particles  converge  along  radial  lines  to 
the  center  of  the  cavity,  compressing  whatever  vapor  might  be  entrapped. 
At  total  collapse,  the  cavity  reaches  a minimum  volume,  and  the  fluid 
particle  motion  is  abruptly  reversed.  A significant  pressure  pulse  is 
then  formed,  and  a second  cavity  develops  which  will  be  smaller  than  the 

v.- 

first.  Higher  pressures  and  extensive  damage  may  well  have  already 
occurred  due  to  projectile  entrance,  exit,  fluid  shock,  and  fluid  pres- 
sures; nevertheless,  the  pressures  associated  with  cavity  collapse  are 
significant. 


80 


AFFDL-TR-75-102 


Williams  (9)  has  measured  the  pressures  at  the  tank  well  due  to 
cavity  collapse.  These  tests  were  for  0.50  cal  AP  at  approximately  0.84 
km/sec  fired  into  three  foot  cubic  tanks.  Kistler  Model  603  accel- 
eration compensated  piezoelectric  pressure  transducers  were  mounted 
flush  with  but  isolated  from  the  tank  wall.  The  loading  on  the  tank 

wall  due  to  the  cavity  collapse  is  shown  in  Figure  24.  These  results 
are  for  the  transducer  at  the  center  of  the  tank  side  wall.  The  other 
transducers,  located  at  the  same  depth  and  halfway  between  the  center 
and  the  edge,  normally  saw  considerably  lower  pressures.  The  t - 0 
reference  for  these  pulses  is  the  time  of  cavity  collapse,  approximately 
26  msec  after  impact.  The  pulses  are  separated  by  about  1.4  to  1.8 
msec. 


There  are  two  items  to  note  from  the  information  in  Figure  24. 
First,  there  was  no  measured  pulse  from  the  cavity  collapse  when  re- 
ticulated foam  was  in  the  tank.  Second,  the  pressure  pulse  was  about  10 
bars  measured  at  the  wall.  It  is  difficult  to  compare  this  pressure  to 
most  recorded  drag  ohase  pressures  since  the  latter  are  normally  taken 
near  the  trajectory.  However,  Williams  also  provides  data  at  the  wall 
for  these  shots.  This  information  shows  pulses  of  shorter  duration  and 
lower  amplitude  and  with  shapes  similar  to  the  pulse  shapes  typically 
measured  near  the  center  of  the  tank  as  shown  in  Figure  25.  The  pulses 
measured  had  peak  pressures  of  about  4 bars  (also  at  the  wall),  and 
pulse  widths  of  about  0.5  msec;  thus  both  in  terms  of  the  peak  pressure 
and  in  terms  of  impulse  per  unit  area,  the  cavity  collapse  pressures  are 
significant. 
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The  analysis  of  the  cavity  expansion  has  been  dealt  with  in  some 
detail  and  it  would  be  desirable  to  use  these  equations  for  cavity 
collapse  also.  However,  Lundstrom's  model  of  the  cavity  is  not  con- 
sidered valid  for  ravity  collapse.  Though  it  has  not  been  verified, 
Yurkovich's  model  is  considered  valid  for  cavity  collapse  provided  care 

is  used  in  applying  the  equation.  Further,  the  cavity  is  observed  to 
become  spherical  before  collapse  begins,  so  that  Yurkovich's  spherical 
bubble  approach  becomes  more  appropriate  here  than  for  cavity  formation. 

5.  DISCUSSION 

Two  principal  relations  were  discussed  in  this  section.  The  first 
was  for  the  transfer  of  energy  from  the  projectile  to  the  fluid,  which 
required  finding  the  terms  in  dEp-dEf+dw.  Here,  the  solution  for 
the  term  dw  presented  no  difficulty. 

The  term  dEp,  however,  was  found  to  depend  upon  knowledge  of  both 
the  projectile  presented  area,Ap(x)  , and  the  projectile  drag  coeffi- 
cient,C0(x)  , where  these  terms  are  given  as  unknown  functions  of  posi- 
tion. The  simplification  is  to  fix  the  presented  area  at  its  value  at 
the  projectile's  maximum  radius.  Then  testing  data  can  be  used  to  infer 
values  to  the  remaining  coefficient  of  drag  function. 

A preferred  approach  would  be  to  observe Ap(x)  , and  then  infer  from 
experiment  the  appropriate  values  of  Cq(6?) , where  8 represents  the 
projectile  attitude,  and  so  also  determines  Ap(0)  . The  tumbling  func- 
tion, Q , would  then  vary  from  test  to  test.  Though  more  difficult  to 
measure,  tumbling  functions  could  then  be  assumed  for  analysis  to  give 
best,  worst,  and  "average"  conditions. 

83 


••  i ' —— amamm r a "n — r r~i  ~ -• 


r'liiifjfc. 


AFFDL-TR-75-102 


The  third  term.dEf  , was  found  by  Yurkovich  and  Lundstrom  through 
very  different  approaches.  Yurkovich  integrated  the  velocity  potential 
for  a spherical  bubble  and  made  the  assumption  that  the  general  form  of 
the  result  would  be  valid  for  similar  shapes.  Lundstrom's  solution 
depends  only  upon  the  assumption  of  an  incompressible  fluid,  and  even 
that  assumption  is  somewhat  mitigated  by  the  finite  upper  limit  of 
integration  in  Equation  131.  Since  Lundstrom's  result  is  more  than 
twice  as  large  as  Yurkovich1 s,  one  of  the  two  should  be  chosen  as  pre- 
ferable. The  Lundstrom  form  is  more  difficult  to  obtain,  but,  on  the 
other  hand  has  been  verified  by  many  tests  (31)  and  has  the  added  ad- 
vantage of  being  derived  from  a physically  more  appealing  model. 

The  second  principal  relation  of  this  section  involves  Bernoulli's 
equation.  Once  solved,  Bernoulli's  equation  gives  the  pressure  field  in 
the  fluid  which  is  used  as  a basis  of  predicting  tank  wall  loading  and 
wall  rupture  due  to  the  fluid  drag  phase  of  hydrodynamic  ram.  In  both 
cases,  Bernoulli's  equation  is  solved  with  recourse  to  a velocity  po- 
tential function.  Lundstrom's  assumption  of  a line  of  sources  leads  to 
a more  complex  model  for  which  no  results  are  claimed  for  cavity  col- 
lapse but  which  has  been  frequently  verified  for  cavity  formation  (31). 
Yurkovich' s model  has  not  been  verified  nor  compared  analytically  to 
Lundstrom's.  But,  it  has  the  advantage  of  being  simpler,  and  has  the 
possibility  of  being  used  for  cavity  collapse,  especially  since  the 
cavity  normally  becomes  spherical  before  it  begins  to  collapse. 
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SECTION  V 

TANK  WALL  LOADING  AND  RESPONSE 

Two  problems  need  to  be  discussed  in  this  section.  Obviously, 
the  response  of  the  tank  wall  to  the  fluid  pressures  needs  to  be  known, 
but  there  are  several  computer  solutions  already  available  which  deal 
with  the  response  of  plates  and  membranes  of  various  shapes  to  given 
loadings.  What  is  needed  more  especially  is  a model  for  predicting 
how  the  fluid  pressures,  which  are  presumed  known  from  the  models  of 
the  previous  section,  load  the  tank  walls. 

Four  approaches  to  loading  or  response  have  been  suggested  by 
authors  working  with  hydrodynamic  ram.  Chou  developed  an  analysis  of 
the  shock  produced  by  hypervelocity  impact.  He  uses  the  IK  yand-Mindlin 
equations  to  describe  the  wall  response  to  fluid  shock.  Yurkovich  did 
not  treat  the  wall  response  in  detail  but  instead  presented  a simpli- 
fied solution  to  illustrate  the  general  considerations  to  be  dealt  with. 
Lundstrom  and  Ball  both  dealt  with  wall  loading.  Lundstrom  presents 
the  Kirchoff  relation  which  should  correctly  predict  wall  loading  but 
which  requires  extensive  computer  time.  Ball  applies  piston  theory  to 
the  fluid  loading  of  the  tank  wall.  Piston  theory  is  presently  the  most 
attractive  approach  because  of  its  simplicity;  none  of  these  solutions 
have  yet  been  verified. 
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1 . FRONT  WALL  RESPONSE  TO  SHOCK 

Chou  (17,  18)  has  referenced  the  analysis  of  Jahsman  (19)  and  ap- 
plied those  techniques  to  the  observations  by  Stepka,  Morse  and  Dengler 
(6)  of  tank  wall  fracture  due  to  impact  by  hypervelocity  particles. 

Chou  presents  the  appropriate  form  of  the  Uflyand-Mindlin  equations, 
gives  the  characteristic  equation  solution,  and  applies  appropriate 
boundary  conditions  for  a circular  hole  in  an  infinite  plate. 


t 

(Assuming  a thin  elastic  plate  with  no  initial  surface  tractions  and 

an  axi symmetric  load,  the  Uflyand-Mindlin  equations  can  be  written  in 

F 

\ polar  coordinates: 

L 
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0R.k2*0fc[*  + |f.] 


(170) 


Most  of  the  terms  in  Equations  168  and  169  are  represented  in  Figure  26 


s'.  I 

i 


Figure  26.  Moments,  Shear  Stress  Resultants, 
Rotations,  and  Displacements  in  a 
Thin  Elastic  Plate 
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The  stress  moments  are  radial  bending, Mr,  tangential  bending, 

M$ , and  twisting,  Mrq.  The  shear  stress  resultants  are  transverse, 
Or  , and  radial ,Qq  ; and  '4'  are  the  rotations  of  the  cross  section 
about  the  tangential  and  radial  axis;w  is  the  transverse  displacement 
of  the  midplane;  and  h is  the  plate  thickness.  Axisymmetric  loading 
conditions  require  that  Mra  "Oj  *^sO  • 

The  following  terms  are  defined  for  Equations  168  and  169.  The 
flexural  rigidity  is  given  by 


0 


Eh3 

12  { I - v*  ) 


(171) 


where  E is  the  modulus  of  elasticity  and  v is  Poisson's  ratio.  A shear 
stress  resultant  constant, k2  , is  introduced  to  account  for  warping  over 
the  plate  cross  section.  Equations  168  and  169  can  be  found  by  inte- 
gration of  Hook's  Law  of  relations. 

Chou  uses  the  method  of  characteristics  to  solve  the  Uflyand- 
Mindlin  equations;  this  method  will  not  be  explained  here,  only  the 
results  will  be  outlined.  One  reference  for  the  method  of  character- 
istics is  Karpp  (20). 
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Substitution  of  Equations  168  and  169  Into  Equations  166  and  167 
yields  a system  of  two  hyperbolic  second-order  equations  inland  w 1 


0 


RR 


L_  * 

c * Vfi 


k2*Gh 


( 4>  + wR)  + -*• 

K R* 


<J> 


(172) 


WRR  “ 


I 

"2  c2 


( 0 + wR) 


-*R- 


F(R,t) 

K2*6h 


(173) 


where  F(R,t)  is  a known  resultant  surface  traction  term,cp  is  the  plate 
velocity,  and  C2  is  the  shear  wave  velocity.  Cp  and  C2  are  defined  by 


E 

/Ml-v*) 


(174) 


C 


c 

2 


(175) 


The  four  physical  characteristics  are  then  given  by 

dR  . r 

W 1 CP 


I± 


(176) 


II 


± . 


dR 

dt 


* ± kg  Cg 


(177) 
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where  the  upper  signs  refer  to  I^and  ligand  the  lower  signs  to  I" 
andll".  Based  on  continuity,  two  final  equations  will  be  written: 


d<J>  = $RdR  + 4>t  dt 


(180) 


dw  ■ wRdR  + dt  (181 ) 

Equations  178  and  179  now  provide  six  equations  in  six  variables: 

4>r  * »*R  » wt  *£>  andw  • 

To  establish  initial  and  boundary  conditions,  recall  that  the  model 
is  an  infinite  plate  with  a circular  hole.  The  radius  of  the  hole  is 
R0  , and  the  plate  is  initially  at  rest  with  no  loads  applied.  Thus, 
initial  conditions  forR0SR<®  at  t = 0 are  £r  * £t  = WR*  wt  = 0 . T 

boundary  conditions  are  £rsmr=0  atR  = R0  . 
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For  additional  comments  and  a discussion  of  the  numerical  tech- 
niques governing  solution  of  this  problem  (especially  the  method  of 
handling  the  discontinuities  in  the  first  derivatives  of  $ andw),  the 
reader  Is  referenced  to  Chou's  works  (17,  18). 

2.  YURKOVICH’ S SPRING  MASS  MODEL 

The  Yurkovich  analysis  of  the  tank  wall  response  (8)  is  intended  to 
provide  insight  rather  than  an  accurate  model,  since  a full  treatment  of 
the  problem  is  both  long  and  complex.  Thus,  with  admitted  simplifi- 
cations, the  problem  can  be  viewed  as  a pressure  loading  on  a circular 
section  of  the  tank  wall.  If  the  fluid  flow  and  the  wall  response  are 
all  in  one  direction,  the  forces  can  be  balanced  relatively  easily  by 
considering  an  equivalent  spring  mass  system. 

With  reference  to  Figure  27,  the  force  equation  can  be  written, 

Ff  • Fw+  F*  (182) 


Figure  27.  Yurkovich' s Wall  Response  Model 
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where  Ff  is  the  force  exerted  on  the  wall  by  the  fluid,  and  Fw  and 
F*  are  the  forces  exerted  by  the  wall;  is  the  force  due  to  the 
kinetic  energy  of  the  wall;  and  Fw  is  the  force  on  the  fluid  due  to  the 
equivalent  spring  displacement.  Before  solving  for  these  forces, 
Yurkovich  assumes  the  wall  is  excited  to  the  first  mode  of  plate  vibra- 
tion, and,  with  reference  to  Timoshenko  (21:452),  approximates  the  wall 
displacement  with 


w • w. 


[•-<£> 


*]' 


(183) 


where  w0  is  the  displacement  of  the  center,  a is  the  radius  of  the  plate 
portion  of  the  wall,  and  R is  the  radial  distance. 

The  force  due  to  the  fluid  is  estimated  with  an  equivalent  force, 
which  is  considered  to  act  at  the  center  of  the  plate;  the  pressure  is 
estimated  by 


P(R,t)  8 Pav  (O 


*o 


where  PQV(t)  is  the  average  pressure  due  to  the  shock, 
force  on  the  wall  by  the  fluid  is 


(184) 
Then  the 


f>  '■[  10  RdRd9'  ■*£  p‘ 


3 'av 


(185) 


where  Equation  188  was  used  for  the  wall  displacements. 
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The  force  on  the  system  due  to  a wall  displacement  of  amount  w 
is  given  by  Fwsk6w  . The  equivalent  spring  constantk®,  is  evaluated 
by  considering  the  potential  energy  of  a compressed  spring,  and  equating 
this  to  the  potential  energy  of  a circular  wall  segment  due  to  bending 
strain  (20:448).  The  potential  energy  of  the  compressed  spring  is  then 


where  Dp  is  the  flexural  rigidity  of  the  plate.  To  evaluate  the  inte- 
gral, Equation  183  is  used  again  for  the  wall  displacement;  after  taking 
the  necessary  derivatives,  expanding,  and  integrating,  the  potential 
energy  becomes 


V = 


7 r 


ZSLY  32 
o / 3 


(187) 


Thus,  for  w=w0  , 


and 
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The  force  due  to  the  kinetic  energy  of  the  wall  can  be  written  in 
terms  of  an  equivalent  mass  of  the  plate,  F*=mew  concentrated  at  the 
center.  The  equivalent  mass  can  be  found  by  evaluating  the  expression 
for  the  kinetic  energy  of  the  plate 

K-  nisi!  s j_  r2irr°  ?hWzRdRt!0  O90) 

K 2 2 Jo  Jo 


where  P , h , and  w are  the  wall  density,  thickness,  and  displacement 
respectively.  The  wall  displacement  is  again  given  by  Equation  183  so 
that  integration  of  Equation  190  then  gives 

«.■*«•(£)  0911 

and 

- ZirP*\az  (192) 

" ' 10 

Expressions  for  all  the  terms  of  the  force  equation  have  now  been 
written:  FO)  , ke  , andme.  Writing  out  the  force  equation  (Equation 

182)  and  rearranging  terms  gives 

«+«>'«=  P°*  — (193) 

•(ft) 
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where 


%* 


094) 


is  defined  to  be  the  natural  frequency  of  the  system. 


A refinement  to  Equation  194  is  possible.  By  considering  all  the 
fluid  within  the  shock  surface  to  not  only  exert  a force  on  the  wall, 
but  to  be  accelerated  to  the  same  velocities  as  the  wall,  an  inertial 


correction  can  be  made  tome.  The  mass  of  the  fluid  within  the  hemi- 

(195) 


spherical  shock  surface  is=^pfa  giving 

3^»h  P*° 


m«  - me  + -j-pf  o -Ztto  [ + -y  J 


where  Pm  is  now  the  density  of  tank  wall  and  p\  is  the  fluid  density. 

Note  that  the  plate  radius, o , must  coincide  with  the  shock  radius, 
r$  , since  more  distant  portions  of  the  tank  wall  are  unaffected  by  the 
shock  pressures. 

The  adjusted  force  equation  now  becomes 


w + (wn)*  w = yTT ( f 1 - — 


(195) 
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and  the  corresponding  natural  frequency  is 


•:3 


(wnV* 


32  De 


(197) 


'•i 


Calculations  for  the  wall  velocities  and  displacements  require  numerical 
integration  of  Equation  196.  The  necessary  difference  equations  are 
derived  by  Yurkovich  (8:41). 

3.  WALL  LOADING  PREDICTED  BY  THE  KIRCHOFF  RELATION 

Neither  Yurkovich' s nor  Lundstrom's  methods  of  calculating  the 
pressures  in  the  fuel  tank  account  for  the  pressures  due  to  reflections 
from  the  tank  walls.  These  pressures  can  be  significant  and  should  be 
accounted  for  before  considering  wall  response,  thus  Lundstrom  presents 
the  Kirchoff  relation  for  calculating  the  pressures  at  the  tank  wall. 

The  pressures  at  each  point  within  the  fluid  as  well  as  at  the  tank  wall 
are  given  by 


P(t,  t)  = P.  (t.t) 


+ 


I 

Air 


I 


+ 
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I 

Air 


/ -L \M 

* ■ \ 

’ • T)  Ids 

Js  r L 

dn  J 

r)  , 1 

dp{t\  r ) i 

1 dr 

— + ~ 

dt  J 

1 dn 

(198) 
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where  r Is  a vector  from  an  arbitrary  origin  to  some  point  where  the 
pressure  is  to  be  calculated* "r  is  measured  from  the  same  origin  to 
successive  points  along  the  tank  wall  from  which  reflected  disturbances 
could  emanate,  and  then  the  source  and  observation  point  are  separated 
by  r = I r — T1 1 , Retarded  time  is  given  by  t , and  n is  the  outward 
normal  to  the  tank  surface.  Equation  197  is  valid  only  for  pressure 

fields  which  satisfy  the  wave' equation,  and  only  for  points  of  the  fluid 
which  are  enclosed  by  surfaces  which  are  concave  inward  (i.e.,  no  part 
of  the  surrounding  surface  can  be  bent  toward  the  point).  Thus,  tank 
walls  with  saddle  points  and  reflections  from  the  cavity  surface  cannot 
be  treated  with  Equation  198. 

dR 

To  apply  Equation  198,  -r—  must  be  written  in  another  form. 

on 

Bernoulli's  equation.  Equation  28,  is  somewhat  simplified  for  no  exter- 
nal forces  and  for  small  velocities  (i.e.,  for  the  point?  not  near  the 
cavity  surface  or  the  bullet): 


d\ 


(199) 


The  partial  derivative  of  Equation  199  with  respect  to  the  outer  normal 
gives 


(200) 
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Lundstrom  anticipates  that  a numerical  solution  of  Equation  202 
will  be  accurate,  but  large  computer  storage  and  long  running  times  have 
discouraged  its  use. 

4.  BALL'S  APPLICATION  OF  PISTON  THEORY 

A final  method  of  solving  the  problem  of  waves  reflected  from  tank 
walls  is  the  use  of  piston  theory  as  advanced  by  Ball  (23).  The  assump- 
tion is  made  that  the  elements  of  the  wall  can  be  treated  as  elements  of 
area  A and  massm  attached  to  springs  of  spring  constant  k as  in  Fig- 
ure 28.  Then  the  total  force  AP(t)  acting  on  the  wall  due  to  the  fluid 
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pressures  is  just  that  which  contributes  to  wall  acceleration  m w,  and 
that  which  contributes  to  the  displacement  of  the  spring  kw  . Thus,  the 
statement  of  the  idealized  one-dimensional  model  is  given  by 


w + kw  = AP  ( t ) 


(203) 


Additionally,  the  incident  pressure  is  modeled  by  a step  pressure 
pulse  beginning  at  time  zero,  as  also  shown  in  Figure  28.  The  pressure 
at  the  wall  is  the  sum  of  the  incident  step  pressure  and  the  reflected 
pressure,  P(t)=Pj  + Pr  (t),  where  the  reflected  pressure  depends  on  the 


ITS"  >"^pri»-.' 
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wall  motion  and  can  be  negative.  The  incident  pressure  is  constant,  and 
the  reflected  pressure  is  approximated  from  the  simplified  form  of 
Bernoulli's  equation,  P(t)=Pj+Pr(t)  , where  the  subscript  r refers  to  the 
reflected  wave.  The  pressure  function  is 


P(»)  * F)  (I) 


(204) 


The  basic  pressure  equation  for  piston  theory  is  just  Equation  204 
with  the  term  <£r(t)  rewritten  in  terms  of  the  fluid  and  wall  velocities. 
To  find  this  form  for<j>r , recall  that^>  satisfies  the  wave  equation. 
Equation  25,  repeated  here  for  convenience: 


c2  v2$  = $ 


(205) 


The  solution  of  the  wave  equation  in  exponential  form  is  given  by 


= aeiM*'  ct)  + /3eiX(*  +c,) 


(206) 


where  a and  f3  are  arbitrary  constants,  and  X is  a scalar  constant  for 
the  wave.  It  follows  from  differentiation  of  Equation  206  that 


u.  = 


du 


(207) 
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and 


ur 


(208) 


The  Incident  and  reflected  pressures  can  now  be  written  from  the  simpli- 
fied form  of  Bernoulli's  equation  and  Equations  207  and  208: 


(j  = A>««, 


(209) 


and 


l>  = -Pcur 


(210) 


Finally,  neglecting  cavitation,  the  wall  velocity  must  equal  the  fluid 
velocity  at  the  wall 


W = U = Uj  + ur 


(211) 


By  combining  Equations  207  and  208  with  the  above  condition,  the  pres- 
sure function  can  be  put  in  the  form 


P ( t ) = + Pc  ( U|  - w) 


(212) 
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The  pressure  function  can  be  written  in  the  final  form  by  consideration 
of  Equations  208  and  209  as 


pressure  in  Equation  203  gives 


where  the  following  terms  were  defined: 


o J 


= 2 f}  - Pew 

(213) 

be  solved.  Using  this 

last  form  of  the 

2 C w + w2  w « y 

(214) 

defined: 

£cA 

2m 

(215) 

II 

+ 

(216) 

2APj 

m 

(217) 

The  solution  of  Equation  214  is  obtained  in  Appendix  D as 


w = 


21_  [(  + -L-  (a"  w,  - a+  «2)] 

- 2£ 


(218) 
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where  the  following  terms  are  defined 


a*  . - i * l 


■ - C - 5 


W,  « • 


a+  t 


w.  = « 


a”  t 


(219) 

(220) 
(221 ) 
(222) 
(223) 


1 1 


I I 


It  follows  that 


-7  (*t  - "8) 


(224) 


and  thus 


Pit)  = 2f»  [ ' - (w,  - H)] 


(225) 


The  only  variable  In  Equation  225  Is  time,  and  most  of  the  parameters 
are  known;  the  spring  constant  is  a parameter  which  must  be  estimated. 
Rather  than  approximate  the  spring  constant  itself,  Ball  estimates 

103 


AFFDL-TR- 75- 102 


the  frequency  of  the  wall  (k/m)  from  the  lowest  natural  frequency  of 
a square  plate  (21:441)  with  width  o and  thickness  h given  by 


i2d-u*y> 


(226) 


where  E Is  Young's  modulus  and  v Is  Poisson's  ratio.  The  incident 

pressure  Pj  is  also  an  unknown  constant,  but  rather  than  stipulate 

P(t) 

values  it  is  more  convenient  to  consider  the  ratio  — — . 


Equation  225  is  plotted  in  Figure  29  for  an  aluminum  wall  2.5  mm 
thick  f or  cu  = 100,  200,  500,  and  1000.  The  problem  divides  naturally 
into  two  parts,  for  early  times  (less  than  10"4  sec)  and  for  later  times 
(greater  than  10”4  sec). 

After  approximating  the  parameters  to  be  encountered  in  hydraulic  ram 
problems,  two  forms  of  the  solution  will  be  found:  the  first  form  will  be 

valid  for  early  times  (less  than  10"4  sec),  the  second  form  will  be  valid 
for  later  times  (greater  than  IQ"4  sec). 


For  this  problem,  a typical  value  of  £ 15  computed  to  be 
7.7  x 104  sec”1  (This  number  is  based  on  densities  of  785  kg/m^  for 

3 

JP-4,  2770  kg/m  for  aluminum,  and  a sound  velocity  of  1360  m/sec  in  the 
fuel.)  From  Equation  226,  the  lowest  natural  frequencies  of  two  half- 
inch plates,  10  and  20  inches  square,  are  found  to  be  600  and  150 

rad/sec,  respectively.  Thus,  the  natural  frequencies  for  this  problem 

* 

can  be  assumed  to  be  less  than  1000  rad/sec  and  £ can  be  expanded  as 
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P = 2770  THICKNESS  = 0.00250 


Figure  29.  Piston  Theory  Prediction  for  Pressure  Ratio, 
P( t)/P-j , at  the  Tank  Wall 


This  expansion  allows  Equation  218  and  Equations  224  and  225  to  be 
rewritten  as: 

v r * “I Tf  « <"2C  +fj£>t  , t99 

w « -4- fl-l-^  + Ue  2^  e 25  1 (zz 
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V 2 1 

For  the  case  where  w /2£  is  less  than  100  sec  , two  time  regimes 
become  apparent  in  Equation  230  and  lead  to  two  approximations: 

For  o < t < I0"4, 


Pltl 

2P, 


(231) 


and  for  |0”4  < t, 


(232) 


These  equations  represent  inertial  and  elastic  effects,  respectively. 

For  w*/2£  = IOO,  the  effects  are  separated  by  about  1 msec. 

Since  the  elastic  response  is  noted  to  occur  relatively  late  in  the 
process,  Ball  next  questions  the  possibility  of  treating  the  fluid 
pressure  pulse  cn  the  wall  as  an  impulsive  load.  For  an  impulsive  load, 
the  incident  pressure  of  infinite  duration  is  replaced  by  an  incident 
pressure  pulse  which  acts  for  timetp  . For  t greater  thantp  , the 
term  P,  in  Equation  217  becomes  zero;  thus,  the  solution  must  satisfy 
two  equations:  For  o<  t < tp, 

w + 2C  w + w*w  - y (233) 
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t 


f t 

f *?.< 


and  for  tp  < t, 


ii  + 2 £ w + (u!w  * 0 


(234) 


Assuming  a solution  of  the  form  w=eat  , and  using  the  same  notation 


as  before, 


a*t  . a't 


w « c,  e + c2  e 


(235) 


w « c,  a • + ct  a e 


_ a"  t 


(236) 


Since  the  functions  w(t)  and  w(1)  are  known  at  t = tp  from  the 
previous  solution,  the  constants  ct  and  c2  can  be  evaluated  by  a simul- 
taneous solution  of  Equations  235  and  236: 


i 1 -a+tp 

T ft  « 


(237) 
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and 

C2  a[w(tp)(-^-1)_^(tp)  a *p  (238) 


Then,  for  £ »a> , the  following  assumptions  are  justified:  £ a£  , 

. «y^  w 2 

a a - jjp- , and  a . With  these  substitutions,  approximate 

values  of  the  constants c , and  c 2 can  be  used  to  write  the  solution  of 
Equation  235,  as  follows: 


(239) 


Finally,  to  describe  an  impulsive  load,  the  finite  pressure  pulse 
is  taken  to  be  large  enough  and  of  short  enough  duration  that  the  wall 
can  be  considered  to  have  remained  stationary  ( w (tp)  = 0 ) while  receiving 
an  initial  velocity  (w(tp)  = /tp)  . It  can  be  seen  from  Equation  217 
that,  y is  a measure  of  the  acceleration  of  the  wall  from  the  incident 
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pressure.  Substitution  of  these  values  for**  and  W into  Equation  239 
gives  the  wall  displacement  due  to  impulsive  loading. 


, „ X'p  2 r 

«*  ( t).  * — • » 

imp  2 t [ 


(240) 


To  determine  the  accuracy  of  the  impulse  load  assumption,  the 
values  of  w(tp)  and  w(tp)  are  obtained  from  Equations  228  and  229  and 

substituted  into  Equation  239.  After  making  these  substitutions  and 

2 

simplifying  the  results  for  the  conditions  thatrl)~ao  and  2)  tp 

2 2£ 

is  small  enough  to  approximate  tp  with  |-£~tp  » the  first  two 
terms  of  the  expansion,  then 


• 1..I  > =,  I's-  ."iT(’"p)  . 

It>tp  2 r 


-2C  v (-2C  +7f)(  1 _tp) 


(241) 


The  first  term  is  a long  time  response,  and  the  second  term  is  a short 
time  response.  Comparison  of  Equation  241  to  240  shows  that  the  impulse 
assumption  correctly  predicts  the  long  time  response,  but  not  the  short 
time  response. 
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5.  DISCUSSION 

Learning  the  response  of  the  tank  wall  is  one  of  the  most  important 
end  products  of  hydrodynamic  ram  studies.  It  is  also  the  most  dif- 
ficult, especially  when  realistic  fuel  tank  geometries  are  of  interest. 
Chou's  approach,  using  the  method  of  characteristics,  seems  particularly 
appropriate  for  the  shock  effects.  This  ir,  because  the  shocks  are  very 
localized,  and  tank  geometries  become  less  important  the  farther  the 
observer  moves  from  the  impact  point. 

All  of  the  other  analyses  were  concerned  with  the  loading  of  the 
wall  to  fluid  pressures.  The  complexities  of  the  problem  become  ap- 
parent in  reviewing  these  analyses.  The  wall  is  loaded  nor.  only  by  the 
pressure  field  created  by  the  slowing  projectile,  but  also  by  reflected 
pressure  from  other  surfaces  as  well  as  the  wall's  own  motion.  Once  the 
loading  on  the  wall  is  known,  this  information  can  be  analyzed  by  sev- 
eral available  wall  response  computer  programs  to  estimate  damage.  The 
Kirchoff  relation  can  be  used  to  compute  the  wall  loading  once  the 
pressure  field  is  known,  but  it  is  a lengthy  process  and,  so,  of  ques- 
tionable practical  value.  Piston  theory  effectively  isolates  the  var- 
iables (Equation  228)  so  that  with  only  the  incident  pressure  field 
known,  the  total  pressure  at  the  wall  can  be  calculated,  and  wall  re- 
sponse can  be  computed. 

There  is  an  added  necessity  for  isolating  the  variables  as  piston 
theory  does.  In  measuring  pressures  due  tc  fluid  drag  on  the  projec- 
tile, the  pressure  field  at  the  wall-fluid  boundary  is  rather  difficult 
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to  obtain.  On  the  other  hand,  the  pressures  in  an  infinite  fluid  (one 
where  the  tank  walls  are  sufficiently  distant  that  there  are  no  reflec- 
tions) are  measured  routinely  and  with  comparative  ease. 

The  flow  of  information  necessary  to  obtain  damage  estimates  is 
indicated  in  a very  general  manner  in  the  following  figure. 


Figure  30.  Information  Flow  for  Predicting  Wall  Response 
Due  to  Hydrodynamic  Ram 
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SECTION  VI 
DAMAGE  REDUCTION 


The  final  goal  for  hydrodynamic  ram  studies  is  to  develop  ways  to 
eliminate  the  extensive  fracture  that  occurs  immediately  after  impact  by 
a projectile.  If  the  damage  could  be  contained  to  a mere  hole  in  the 
entrance  and  exit  walls,  that  level  of  damage  would  be  considered  ac- 
ceptable and  the  fuel  loss  would  be  minimal.  Four  methods  to  reduce 
hydrodynamic  ram  damage  have  been  tested,  but  only  one  method,  foam- 
backed  walls,  has  shown  much  promise.  Additionally,  there  are  three 
analysis  techniques  that  can  be  applied  to  the  latter. 


1.  DAMAGE  REDUCTION  TECHNIQUES 

Two  reports  treat  damage  reduction  methods  in  detail;  the  first,  by 
Stepka  and  Morse  (25),  deals  with  ways  to  minimize  damage  to  spacecraft 
fuel  tanks  through  fuel  pressure  attenuation.  Bristow  and  Lundberg  (26) 
consider  the  specific  problem  of  bullets  impacting  aircraft  fuel  tanks. 
Damage  reduction  is  also  discussed  in  several  other  reports  and  has  been 
the  subject  of,  or  included  in,  several  testing  programs. 


In  most  damage  reduction  techniques,  the  goal  is  to  absorb  the 
fluid  energy  before  it  can  damage  the  tank  walls  - reasoning  that  to 
minimize  damage,  the  shock  or  high  pressure  pulse  must  be  spread  out 
over  time  or  space  or  both.  The  suggested  methods  include  adding  sub 
stances  to  the  fuel  to  retard  the  propagation  of  the  wave,  and  adding 
shock  absorbing  material  to  the  inside  of  the  tank  wall. 
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An  entirely  different  approach  was  taken  by  Bristow  and  Lundberg 
(26).  They  tested  a viscoelastic  additive  to  the  fuel  whose  viscosity 
increases  with  increasing  strain  rate.  The  high  pressure  regions  gen- 
erated by  the  projectile  were  expected  to  generate  a highly  viscous 
region  and  thus  increase  the  drag  on  the  projectile. 

To  test  this  concept,  0.50-cal  projectiles  were  fired  into  a tank 
containing  a mixture  of  fuel  and  up  to  1%  viscoelastic  additive.  Un- 
fortunately, there  was  no  measurable  change  in  the  hydrodynamic  ram  pa- 
rameters; projectile  characteristics,  fluid  drag,  and  pressure  pulses 
remained  the  same.  Thus,  the  addition  of  a viscoelastic  substance  to 
the  fluid  in  this  manner  has  not  affected  the  hydrodynamic  ram  damage 
mechanism. 

As  an  addition  to  this  testing  program,  Bristow  and  Lundberg  added 
ordinary  whiffle  balls  to  the  fuel  to  see  if  they  had  any  retarding  effect 
on  the  pressure  pulse.  But  again,  there  was  no  measurable  change  in  the 
hydrodynamic  ram  characteristics. 

A second  ram  attenuation  technique  is  to  absorb  the  fluid  energy 
before  it  impacts  the  tank  wall  by  the  addition  of  a gas  to  the  fuel. 
Williams  (9)  recommended  that  nitrogen  be  bubbled  through  the  fuel. 

Clark  (27)  later  investigated  the  attenuation  properties  of  air  en- 
trapped in  water  by  foam  and  found  significant  pressure  reduction;  in 
addition  he  noted  that  the  ram  attenuation  properties  of  freon  should 

be  superior  to  air.  There  are  two  advantages  to  using  the  gas: 

(1)  the  bubbles  are  not  needed  until  the  aircraft  is  near  the  target, 
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so  that  the  generators  need  not  be  turned  on  for  most  of  the  mission  and 
need  not  even  be  provided  for  all  the  tanks,  (2)  inerting  of  the  ullage 
is  an  advantageous  byproduct;  both  freon  and  nitrogen  are  inert  and, 
further,  nitrogen  will  scrub  the  oxygen  from  the  fuel.  Williams  has 
recommended  that  the  volume  of  the  gas  be  about  6%  N,  which  is  difficult 
to  attain  with  available  generators. 

The  most  promising  method  to  reduce  the  fracture  at  the  entrance  and 
exit  walls  is  to  back  the  tank  wall  panels  with  a foam  which  isolates  the 
fluid  from  the  wall.  Stepka  and  Morse  (25)  studied  this  technique  for 
application  to  spacecraft  fuel  tanks,  and  more  recently  this  same  con- 
cept was  tested  for  use  in  aircraft  fuel  tanks  (28). 

A typical  foam-backed  panel  is  just  an  ordinary  tank  wall  backed 
with  an  easily  compressible  material  such  as  polyurethane  foam.  Typical 
configurations  for  test  purposes  are  illustrated  in  Figure  31.  The 
plastic  sheet  between  the  foam  and  the  fluid  allows  the  easily  compres- 
sible material  to  isolate  the  tank  wall  to  some  extent  from  the  fluid 
shock  pulses.  Several  materials  besides  the  polyurethane  foam  have 
been  tested,  including  air,  rigid  foams,  and  honeycomb  materials. 

Stepka  and  Morse  conducted  several  tests  using  different 
hypervelocity  projectiles  against  several  target  configurations.  Appro- 
priate parts  of  those  tests  have  been  extracted  and  reproduced  in 
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Figure  31.  Test  Configurations  for  Ram  Attenuating  Foams 


Table  I.  These  results  indicate  that  both  foam  and  honeycomb  materials 
provide  significant  protection  of  the  fuel  tank  wall  from  hypervelocity 
impact.  This  moderate  success  implies  that  some  method  of  isolating  the 
fuel  from  the  tank  wall  can  at  least  reduce  the  degree  of  fracturing, 
even  for  ordnance  impact  into  fuel  tanks. 


Dynamic  Science  (28)  conducted  tests  of  the  same  attenuation  con- 
cept using  14,5  mm  AP1-T  against  a representative  integral  fuel  tank. 
Information  was  obtained  from  two  vertical  panels  and  from  the  bottom 
panel  for  a test  configuration  similar  to  B of  Figure  31,  except  that 
the  panels  were  reinforced  with  Z-stiffeners.  Both  flexible  and  rigid 
foams  were  tested.  A cursory  examination  of  the  damage  data  from  this 
program,  summarized  in  Table  II,  indicates  that  a shock  absorbing  mate- 
rial can  significantly  reduce  the  damage  to  the  front  wall,  but  that 
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TABLE  I 

FRACTURE  AND  PUNCTURE  DATA  FOR  SOME  RAM  ATTENUATING 
FOAMS  FROM  HYPERVELOCITY  IMPACT 

B1  - air  Cl  - 3/4  in.  open  or  closed 

B2  - fibrous  matting  cell  polyurethane  foam 

B3  - foam  rubber  D2  - 1 in.  honeycomb 

B4  - 1/2  in.  open  cell 
polyurethane  foam 
B5  - 3/8  or  3/4  in.  closed 
cell  polyurethane  foam 
B6  - 1 in.  honeycomb 


Projectile: 

Nylon  cylinder  (5.56  mm  length  and  diameter) 

configuration 

velocity  (km/ sec) 

fracture  (f)  or  puncture  (p) 

A1 

4.95,  6.01,  6.43,  6.46 

f,  f,  f> 

f 

A2 

5.48,  5.73 

f,  f 

B1 

5.94 

f 

B2 

5.88,  6.14 

f 

B3 

6.28 

P 

B4 

5.63,  5.69,  5.82,  6.09 

P,  f,  p» 

p with  small  cracks 
and/or  deformation 

B5  (3/8) 

6.00,  6.09 

f,  f 

B5  (3/4) 

5.63,  5.78,  5.87,  5.97 

P.  P»  P» 

P 

B6 

5.C7,  5.94 

P»  P 

Cl 

5.94,  5.97 

P,  P 

C2 

5.48,  5.82,  5.94 

P,  P,  f 

116 


A1  - basic 
A2  - double  wall 
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TABLE  11 

DAMAGE  SUMMARY  FOR  A TANK  IMPACTED  BY  14.5  API-T 


Projectile:  Soviet  14.5  mti  API-T,  velocity  range  997  to  1050  m/sec 

Tank  Dimensions.  3 ft  x 3 ft  with  either  a 1 ft  or  2 ft  shotllne 

Test  Panels:  18  In.  x 18  In.  with  4 Z-stlffeners,  each  stiffener  had  18  rivets 


1 

ft.  shotllne 

2 ft  shotline 

panel  damage'  * 

panel  damage1 

Impact  vel . 

(rivets  rema1n1ng)‘ 

Impact  vel. 

(rivets  remaining) 

1 

(km/ sec) 

entrance 

exit 

(km/sec) 

entrance 

exit 

bottom 

A basell 

ne 

1030 

mod(26)» 

des(13)* 

1043 

sev (27)* 

sev(15)* 

(67)* 

81 

1.5  In.  fle> 

Ible  foam 

1024 

m1n( 1 5) * 

des(O)* 

1034 

min(17)* 

sev(8)* 

(34)* 

1027 

ml  n ( 1 9 ) * 

des(O)* 

- 

m1n(14)* 

sev(16)* 

8 * 

1029 

m1n(30)‘ 

des(8)* 

“ 

sev(19) 

sev(19)* 

(10)* 

82  1.5  In. 

rigid  foam  ove 

r 1.5  In  flexible  foam 

1017 

m1n(37) 

sev(12)* 

- 

m1n(28)* 

sev(17)* 

(1°)(2) 

997 

min  40) 

sev(  6)* 

1028 

min  19)* 

ext  20)* 

23)0) 

1026 

min(13)* 

sev(l3)* 

min(l8)* 

ext (22)* 

(23) 

B3  3 in.  flex 

Ible  foam 

1050 

m1n(3l )* 

sev(13)* 

1043 

m1n(23)* 

mod(56) 

(101* 

1041 

min(17) 

des  10  * 

- 

min(18)* 

mod  10 

21  * 

1038 

min(25) (2) 

sev(ll  )* 

1041 

min(21  )* 

mod(20)* 

(19) 

B4  3 in.  rl 

gld  foam 

1023 

mod(40) 

sev(14)* 

1026 

m1n(32)* 

ext (45) 

(401 

1024 

m1n(30) 

des ( 6)* 

1023 

min  23  * 

sev(2l)* 

(70) 

1026 

m1n(55)* 

sev (20)* 

” 

m1n(29)* 

mod(20)* 

(71) 

a 

3 


4 :i 

V?  * 


gfl 


min  « minor  damage  - projectile  silhouette 
mod  ■ moderate  damage  - cracks 

ext  * extensive  damage  - large  cracks,  some  rupture,  80%  of  panel  gone 
sev  * severe  damage  - large  ruptured  area 
des  1 destroyed  panel 

* « all  revets  failed 

(N)  5 N Is  the  number  of  rivets  that  remained  In  the  panel 
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this  same  material  may  have  no  marked  effect  on  protecting  the  rear 
wall.  It  must  be  added  that,  most  frequently,  combat  inflicted  damage 
is  from  an  impact  into  the  bottom  of  the  tank,  and  that  a flexible  foam 
then  will  already  be  compressed  to  some  extent  due  to  the  weight  of  the 
fuel,  thus  reducing  its  ability  to  absorb  shocks. 

2.  DAMAGE  REDUCTION  ANALYSIS 

There  are  three  analyses  of  compressible  foams  available.  Stepka 
and  Morse  (25)  used  an  analysis  based  on  shock  Hugoniot  curves  to  ap- 
proximate the  maximum  pressures  at  the  interfaces  of  the  different  ma- 
terials, but  that  analysis  does  not  account  for  the  thickness  of  the 
materials.  Torvik  (29)  has  suggested  an  analysis  based  on  the  "snow- 
plow" effect  which  will  approximate  the  reduction  of  the  pressure  as 
the  pulse  traverses  a distance  of  foam.  Williams  (9)  gives  a calcula- 
tion for  the  minimum  thickness  of  foam  based  on  the  time  required  for 
the  rarefaction  wave  to  overtake  the  shock  front. 

a.  Shock  Hugoniot  Analysis.  A minimum  amount  of  information 
about  shock  Hugoniots  is  presented  here  to  aid  in  using  this  analysis 
method.  For  further  information  the  reader  is  referred  to  the  litera- 
ture, especially  Duvall  (30)  and  Cole  (3). 

For  this  shock  analysis,  the  effects  of  the  projectile  and  the 
hole  it  creates  are  ignored.  The  shock  pulse  is  assumed  to  be  created 
in  the  fluid  through  a nearly  instantaneous  deposition  of  energy  by  the 
projectile.  This  analysis  considers  only  the  maximum  pressures  impacting 
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at  the  fluid-foam  boundary  as  described  by  a one-dimensional  model. 
Boundary  conditions  determine  the  state  of  the  shock  which  can  exist 
at  the  boundary,  and  hence  the  shock  which  can  initially  be  transmitted 
through  the  foam.  The  shock  can  be  described  similarly  at  the  foam-tank 
wall  boundary. 

The  set  of  states  that  can  exist  at  the  material  interfaces  is  de- 
scribed by  shock  Hugonlot  curves.  These  curves  represent  the  locus  of 

n 

all  states  of  a material  that  can  be  reached  from  given  initial  condi- 
tions. In  a fluid,  the  Hugoniot  is  a plot  of  the  pressures  against 
particle  velocities;  in  a metal,  the  pressure  is  replaced  by  the  stress 
normal  to  the  shock  front.  The  information  necessary  to  generate  the 
Hugoniot  must  be  obtained  from  experiment,  but  at  the  same  time  this 
data  allows  the  development  of  an  equation  of  state  from  which  the 
Hugoniot  can  later  be  extracted. 

When  a shock  wave  reaches  a boundary,  part  of  it  is  transmitted, 
and  part  is  reflected.  The  transmitted  wave  will  also  be  a shock,  but 
the  reflected  wave  may  not  be.  The  reflected  wave  may  be  a compressive 
wave  or  a rarefaction  wave  depending  on  whether  its  pressure  is  higher 
or  lower  than  the  incident  wave.  If  it  has  a higher  pressure,  it  is 
compressive,  and  may  be  a shock.  The  Hugoniot  curve  in  Figure  32  may 
then  be  reflected  about  the  vertical  axis  and  its  intersection  with  the 
P * 0 axis  aJjusted  to  describe  the  states  available  to  the  reflected 
wave,  I . If  the  reflected  wave  is  a rarefaction  wave,  the  reflection 
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Figure  32.  Direct  Hugoniot,  Reflected  Hugoniot,  and 
Reflection  Adiabat 


adiabat  must  be  used  but,  in  practice,  this  curve  is  nearly  identical  to 
the  continuation  of  the  reflected  Hugoniot,  and  so  the  reflected  Hugoniot 
is  used  exclusively. 

The  intersection  of  the  reflected  Hugoniot  with  the  particle  ve- 
locity axis  is  the  particle  velocity  Up.  For  the  shock  wave  transmitted 
to  the  fluid  from  the  projectile.  Up  is  taken  to  be  the  velocity  of  the 
projectile  immediately  after  it  has  penetrated  the  tank  wall.  For  thin 
walls,  the  velocity  may  be  taken  as  the  impact  velocity.  After  locating 
the  intersection  point  on  the  particle  velocity  axis,  the  reflected  and 
direct  Hugoniots  will  intersect  at  some  states,  which  describes  the  only 
values  of  fluid  pressure  and  fluid  particle  velocity  the  curves  can 
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share.  This  is  illustrated  in  Figure  33.  The  Hugoniots  of  additional 
materials  can  then  be  added,  and  a final  estimate  of  the  expected  maxi- 
mum pressures  can  be  made. 

Figure  34  illustrates  how  the  maximum  expected  pressures  at  the 
fuel  tank  wall  are  determined  graphically  for  layers  of  different  mate- 
rials. 

b.  “Snowplow"  Analysis.  Following  the  works  of  Hermann  (31)  and 
Salvadori,  et  al  (32),  Torvik  (29)  has  described  the  application  of 
"snowplow"  model  to  the  attenuation  of  shock  waves  in  porous  materials. 
This  model  predicts  the  decay  of  the  peak  pressure  with  the  square  of 
the  distance  the  shock  wave  has  penetrated. 


Figure  33.  Identification  of  Boundary  and  Initial  Conditions 
Using  the  Direct  and  Reflected  Hugoniot 
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Torvik  considers  a porous  material  (foam)  in  one  dimension  as  it  Is 
being  subjected  to  a shock  wave,  as  illustrated  in  Figure  35.  The  foam 
has  been  approximated  as  a "locking  solid,"  i.e.,  any  applied  pressure 
compresses  the  foam  its  maximum  amount.  Thus,  before  the  shock  front, 
the  foam  has  a density  of/50,  and  immediately  behind  the  shock  front, 
the  foam  has  a density  of  . The  model  further  requires  constant 
momentum  in  the  material  behind  the  shock  front.  Thus,  as  the  mass  of 
the  foam  is  increased  with  the  passing  shock  front,  the  particle  ve- 
locities must  decrease  accordingly.  At  the  same  time,  the  particle 
velocities  are  required  to  be  uniform  so  that  p f is  constant  throughout 
the  compressed  region. 


SHOCK 

FR0!«T 


Figure  35.  "Snowplow"  Diagram 
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The  following  terms  can  be  defined  with  reference  to  Figure  35: 
u,  is  the  particle  velocity  just  behind  the  shock  front;  then,  with 
the  assumption  of  a locking  solid,  u,,  the  particle  velocity  for  all  the 
material  behind  the  shock  front,  is  just  u,.  The  shock  moves  with 
velocity,  Us,  and  has  traversed  a distance  I into  the  foam.  An  arbi- 
trary particle  of  the  compressed  foam  is  located  at  x. 


The  pressure  is  to  be  evaluated  ultimately  from  an  integration  of 
the  conservation  of  momentum  (Equation  2)  written  in  two  dimensions: 


dp  _ du 

dx  ~ ~ po  dt 


(242) 


The  term is  just  4^-since  the  locking  solid  assumption 
requires  that  all  the  particle  velocities  be  the  same,  and is  the 

du,  0 

velocity  of  the  shock  front,  U8.  To  evaluate-^,  U,  is  determined  from 
the  impulse  per  unit  area  I0.  Using  £ as  a dummy  variable  for  inte- 
gration, a slice  of  the  foam  has  mass  per  unit  area/&0d£.  Thus, 


= u,  IX) 


(243) 


is  integrated  easily  to  give 


I0  * 


Po 


(244) 
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The  desired  partial  derivative  is  then  given  by 


The  pressure  follows  from 


(245) 


* -r%-  u.  "* 


(246) 


Finally,  the  conservation  of  mass  across  the  shock  front  (Equation  12) 
can  be  rearranged  to  give 


_ _ Pi~P0 

ul  * “ — 

Pi 


(247) 


where  u,  was  defined  earlier  (Equation  31)  as 


'■  “ U, 


(248) 


Taking  u,  as  constant  for  the  locking  material,  and  using  Equation 
242  for  u,  allows  the  integration  of  Equation  246  resulting  in 


P = — =r 


(249) 
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The  assumption  of  the  locking  solid  leads  to  the  prediction  of  a 
triangular  shock  pulse  shape  shown  in  Figure  36,  with  a maximum  ampli- 
tude determined  by  the  distance  the  shock  has  traveled  in  the  fu.m. 

c.  Annihilation  Calculation.  Another  method  of  analysis  can  be 
applied  to  the  foam.  Williams  (9:16)  observed  that  after  the  shock  wave 
has  progressed  through  the  foam  a short  distance,  a rarefaction  wave  is 
formed  behind  the  shock.  The  rarefaction  wave  has  a higher  velocity 
than  the  shock  wave  and  so  overtakes  it;  since  the  particles  behind  a 
strong  rarefaction  wave  have  very  low  densities,  the  shock  can  be 
annihilated. 

The  shock  is  assumed  to  be  a pulse  of  duration,  tp  . During  this 
time  the  shock  progresses  the  distance  uttp  , while  the  free  surface 
travels  only  utp.  The  rarefaction  wave  is  then  formed,  and  will  catch 
the  shock  wave  in  the  time  tc  . The  rarefaction  wave  travels  with  ve- 
locity Cr,  the  release  velocity.  The  relationship  between  these  dis- 
tances is  shown  in  Figure  37.  Conservation  of  mass  across  the  shock 
(Equation  12)  is  used  to  write  the  particle  velocity,  u,  in  terms  of  the 
shock  velocity 

(250) 


u = u.  ( I 


' V") 


Using  the  distances  shown  in  Figure  31  the  following  relationship 
can  be  made: 


u 


t 


+ Cr  tc 


Mil  I Tl 'iiii'i ThfiaKii'i’  if  "Mr 


(251) 
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v-  Thus,  the  distance  required  for  annihilation  of  the  shock  is 
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3.  DISCUSSION 

It  was  noted  In  this  section  that  Williams'  description  of  shock 

annihilation  predicts  the  total  elimination  of  the  shock  within  only  0.3 

cm  of  foam.  Torvik's  snowplow  model  also  can  be  used  to  estimate  the 

distance  required  for  a stipulated  amount  of  pressure  reduction,  and  a 

comparison  of  the  two  predictions  is  useful.  A rough  number  value  for 

the  snowplow  prediction  can  be  given  by  using  Equation  249  with 

x « X (the  maximum  pressure  plane).  The  same  variables  as  In  the 

3 3 

Williams'  calculation  are  used:  p0  = 0.186  g/cm  , p^  ~ 0.833  g/cm  , 

+p  - 8 ysec,  and  an  initial  pressure  of  400  bar  (5800  psi)  - approxi- 
mately the  peak  measured  pressure  1/2  inch  from  the  impact  point  of  a 
0.50  cal  projectile.  Thus  2.66  cm  are  required  to  reduce  the  pulse  to 
10  bars,  and  8.42  cm  are  required  to  reduce  the  pulse  to  only  1 bar. 


!'JJUi,'l.l.llU>  lljUpjlll'LWp! 


The  snowplow  model  may  accurately  predict  the  pressure  reduction 
due  to  the  impulsively  formed  shock,  but  it  does  not  predict  any  change 
in  the  impulse  delivered  to  the  wall  for  typical  fuel  tank  ram  problems. 
Consider  the  following  simplified  example.  Let  the  400  bar,  8 ysec  im- 
pulse above  be  reduced  to  one  bar  with  about  9 cm  of  foam.  If  the  im- 
pulse remains  unchanged,  then  the  one  bar  pulse  is  spread  out  to 
_2 

0.32  x 10  sec.  But  from  the  section  on  wall  loading,  the  natural 

_2 

frequency  of  the  wall  was  found  to  be  of  the  order  of  one  cycle  per  10 
sec.  Thus  the  wall  still  sees  the  load  as  nearly  impulsive.  Since  the 
wall  sees  both  loads  as  nearly  identical  impulses,  the  wall  response  is 
not  much  affected  by  the  reduced  pressure. 
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Tests  to  date  of  up  to  3 inches  of  foam  have  shown  only  moderate 
damage  reduction,  and  that  at  the  front  wall  where  the  shock  would  also 
be  expected  to  predominate.  This  could  indicate  that  the  shock  Is 
predominant  in  the  early  reaction,  and  that  the  foam  has  some  moderating 
effect  on  the  shock.  This  test,  though,  does  certainly  Indicate  the 
importance  of  testing  foams  and  foam  configurations  for  their  responses 
to  shocks  and  pressures  separately. 


Nitrogen  is  currently  being  tested  as  a method  to  Inert  the  ullage. 
Nitrogen  gas  bubbles  in  the  fuel  would  not  only  scrub  the  oxygen  but 
would  also  provide  ram  attenuation.  Studies  of  this  method  must  also  be 
recommended . 


Finally,  a remark  about  self-sealing  materials  is  appropriate. 
Although  primarily  intended  to  seal  the  hole  after  the  projectile  has 
passed  through,  these  materials  also  exhibit  ram  attenuation  properties. 
Self-sealing  materials,  when  combined  with  other  attenuation  materials, 
may  be  particularly  desirable.  The  preceding  analyses  should  be  appli- 
cable after  appropriate  modification. 
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SECTION  VII 
CONCLUSIONS 


The  ultimate  objective  of  hydrodynamic  ram  studies  is  to  elimi- 
nate the  extensive  fracture  and  rupture  that  follows  the  Impact  of  a 
fuel  tank  by  an  ordnance  projectile.  Although  this  goal  may  eventually 
be  shown  to  be  unreachable,  it  is  certainly  possible  to  reduce  damage 
to  a minimum.  To  do  this,  it  Is  necessary  that  each  phase  of  the 
phenomenon  be  observed  in  detail  and  fully  analyzed.  An  understanding 
of  the  phenomenon  will  permit  the  accurate  evaluation  of  ram  attenuat- 
ing mechanisms  and  may  result  In  suggestions  for  others.  More  impor- 
tantly, this  understanding  of  the  ram  event  will  allow  the  analysis  of 
fuel  tanks  and  fuel  tank  systems  while  they  are  still  in  the  planning 
stages. 


Thus  far,  analysis  and  testing  of  liquid-filled  containers  have 
shown  that  hydrodynamic  ram  is  a very  complex  phenomenon,  and  have 
tried  to  simplify  its  study  by  separating  it  into  its  major  events. 

This  simplification  will  prove  most  useful  if  fracture  can  be  associ- 
ated primarily  with  only  one  of  these  major  events.  It  remains  now  to 
be  demonstrated  when  and  how  fracture  occurs  so  that  it  can  be  deter- 
mined which  phase  predominates  in  causing  rupture  under  given  conditions. 

Some  aspects  of  hydrodynamic  ram  were  treated  rather  fully. 

Yurkovich  gave  an  analysis  of  the  shock  phase  and  the  drag  phase,  and 
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an  approximate  approach  to  the  wall  response.  The  shock  phase  work 
has  been  verified  only  against  data  from  hypervelocity  studies  and  its 
applicability  to  ordnance  projectiles,  while  reasonable  as  a limiting 
case,  has  not  been  verified. 

Lundstrom  has  described  the  drag  phase  of  the  event  in  considerable 
detail  and  developed  a solution  of  Bernoulli's  equation  to  describe  the 
pressure  field. 

Chow  has  used  the  method  of  characteristics  to  describe  tank  wall 
response  to  shocks,  but  his  work  has  not  been  verified  for  bullets. 

Ball  has  recommended  piston  theory  for  studying  the  wall  response  to 
both  fluid  shocks  and  pressures. 

The  next  logical  step  in  the  hydrodynamic  ram  studies  is  to  veri- 
fy existing  models  and  their  assumptions.  Shock  fronts  at  the  impact 
point  and  pressure  pulses  in  the  fluid  have  both  been  measured  routinely. 
This  information  has  permitted  the  calculation  of  wall  loading,  which 
is  necessary  to  predict  wall  response.  But  because  the  pressure  load- 
ing at  the  wall  has  been  difficult  to  measure,  the  wall  loading. models 
have  not  been  verified.  However,  Ytterbium  gages  can  be  used  to  get 
pressure  measurements  at  the  wall-fluid  interface,  though  with  some 
difficulty.  In  addition,  wall  response  itself  has  not  been  measured 
except  at  isolated  points  with  strain  gages;  yet  wall  response,  includ- 
ing fracture,  can  be  measured  over  areas  larger  than  2 sq  ft  and  at 
time  intervals  of  approximately  lOy  sec  using  moire  fringes  and  high 
speed  photography. 
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To  a large  extent,  this  lack  of  data  exists  because  It  Is  very 
difficult  to  obtain.  Ytterbium  gages  have  been  used  to  get  pressure 
measurements  at  the  wall  fluid  Interface,  but  only  with  much  diffi- 
culty. A moire  fringe  technique  to  measure  wall  motion  has  been  sug- 
gested recently,  but  the  system  is  critically  dependent  upon  finding 
an  adequate  light  source. 

Although  many  pieces  of  information,  such  as  strain  and  strain 
rates  within  the  tank  walls,  may  be  essential  to  predicting  the  response 
of  the  tank  to  bullet  impact,  they  cannot  be  measured  directly.  This 
information  must  be  properly  inferred  from  verified  analytical  models 
or  from  computer  codes  whose  accuracy  is  known.  Thus  it  is  expected 
that  the  abundance  of  detail  available  from  two-dimensional  wave 
propagation  codes  will  make  those  codes  invaluable  for  understanding 
the  damage  mechanisms,  failure  criteria,  and  attenuation  possibilities 
associated  with  hydrodynamic  ram. 

The  study  and  prediction  of  wall  response  is  of  fundamental  im- 
portance. Once  the  loading  is  determined,  several  codes  can  be  used 
to  predict  deformation  and  failure.  STAGS  and  BR-1  have  both  been 
modified  to  account  for  hydrodynamic  ram  loading.  But  accurate  load- 
ing predictions,  particularly  at  the  exit  and  entrance  regions,  are 
not  available  from  observation  or  current  models.  Only  the  general 
Lagrangian  or  Eulerian  codes  mentioned  above  are  able  to  predict  this 
type  of  loading. 
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The  following  tasks  and  study  areas  are  recommended. 

1.  Collection  of  wall  fracture  data  to  determine  the  extent  of 
wall  fracture  as  a function  of  time.  This  data  should  be  obtained  for 
several  projectiles,  velocities  and  masses,  and  for  various  degrees  of 
tumble  at  entry.  This  data  needs  to  be  obtained  for  both  entrance  and 
exit  walls;  similar  data  needs  to  be  obtained  for  side  walls. 

2.  A model  needs  to  be  developed  for  exit  wall  impact. 

3.  The  existing  fluid  shock  model  must  be  verified  for  ordnance 
projectiles. 

As  a final  comment,  the  explication  of  available  computer  pro- 
grams in  study  areas  of  hydrodynamic  ram  has  been  beyond  the  intent  of 
this  report.  Yet,  the  use  of  these  codes  is  essential;  nearly  all  of 
the  preceding  analysis  has  been  programmed  to  facilitate  obtaining 
numerical  results.  Additionally,  most  of  the  work  is  so  complex  that 
it  can  only  be  approached  with  the  aid  of  the  computer. 
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APPENDIX  A 

GRADIENT,  DIVERGENCE,  CURL,  AND  La P LAC I AN  EQUATIONS 
IN  CYLINDRICAL  AND  SPHERICAL  COORDINATES 

cylindrical  , spherical 
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APPENDIX  B 

EXAMPLES  SHOWING  THAT  THE  SCALAR  POTENTIAL  SATISFIES 
THE  WAVE  EQUATION 


To  show  that  the  scalar  potential  <f>  satisfies  the  wave  equation, 
Equation  1 is  rewritten  as: 


7 


I *P 


It  is  necessary  to  first  rewrite  the  differential -g^-.  Density  is 
a function  of  both  pressure  and  entropy,  but  entropy  can  be  assumed 
constant:  this  assumption  is  verified  numerically  by  Cole  (3:sec.  2.6). 

Using  the  definition  of  sound  velocity  c*  « and  assuming  constant 

dP  dP  dp  ^ ® 

entropy,  the  derivative  dt  8 ~ "gf  leads  to 


The  assumption  of  constant  entropy  also  allows  (from  di*TdS  + TdP ) 


7 p = p 7i 


and  thus  the  conservation  of  momentum  equation  (Equation  2)  becomes 


, _vi 

dt 
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When  written  out  and  simplified,  this  last  equation  is 


du 

dt 


- Vi  +tx  (Vxtf) 


(B5) 


where  Vi  s7l-(u-7)  u-uX(Vxu)  follows  from  the  definition  of 

■A 

kinetic  enthalpy:  I * i -i0+  . The  assumption  of  irrotational 

motion  then  eliminates  the  last  term  on  the  right  and  gives  ~-V<£*-Vl 
which  leads  to  the  conclusion  that 


I * - 


d<p 

dt 


(B6) 


Now  the  expression  for  V ♦ *u  can  be  expanded  to  give 


---if  Hi~+ “ v') 


(B7) 


where  the  two  terms  in  parenthesis  are  to  be  rewritten.  From  the 
kinetic  enthalpy: 


_£L_  s 

at 


(B8) 


= -jf-  - i *7T  ( T “ ) + i(o  V)(if  0 ) 
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APPENDIX  C 

DERIVATION  OF  THE  CONSERVATION  EQUATIONS  IN  NONDIMENSIONAL  FORM 
Derivation  of  the  Conservation  equations  in  nondimensional  form. 
1.  CONTINUITY  EQUATION 

The  continuity  equation  is  given  by  Equation  6 as 


dt  r dr  dr  r 


(Cl) 


Using  the  substitutions 


P SP0  P (7  , M8) 


S Us  u ( r , Ms  ),  and  r * r,  7 (C2) 


leads  to 


(JL  JL.  ♦ +p  u~  + 

\ dr  dt  dM*  dt  / 0 \ dr  dr  dMf  dr  / 

(C3) 

► „ u.,  JL  + JL  ) + ifi££JL.0 

r°  8 \ dp  dr  dMt  dr  / r%  r 
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APPENDIX  D 

THE  SOLUTION  OF  THE  FORCE  EQUATION 


The  solution  of  the  force  equation 

w+2£w+u j2  w = y (Dl) 


is  a little  lengthy,  but  straightforward.  It  is  a second  order  linear 
differential  equation  that  can  be  handled  by  the  method  of  variation  of 
parameters,  as  in  Ross  (24).  Assuming  a solution  of  the  homogeneous 
part  of  the  form  w«  ea*  yields  a complementary  solution 


*c 


c,  • 


a+t 


a t 


(02) 


where 


<•*  * - C + £ * (i») 

a"  * - C - C * (04) 

C *•  (05) 


Then,  by  the  variation  of  the  parameters  method,  the  following  equation 
remains  to  be  solved: 


144 


W|  = e 


<*+t 


m9  a • 


a t 


The  parameters  V|  , and  V2  are  to  be  determined  from 


(D7) 


(D8) 


(D9) 


and 


_ f /"  i 
v2  ~ j W 


dt 


(DIO) 


where 


W * 


»«  *2 


= 2«|"2C 


(Dll) 


The  solution  of  the  particular  equation  now  yields 


v\  *i  + v2w2 s — r 

to» 
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The  constants  c,  and  ct  of  the  complete  solution  are  determined  from  the 
Initial  conditions,  w«  w»0  at  t»0.  Thus, 


a 

•I 


ya' 


2w*5* 


(013) 


and 


r<* 

c.  J T 

a 2w*C* 


(014) 


I 


Making  these  substitutions  in  Equation  D6  gives  the  solution  of  the 
force  equation 


w » 


One  differentiation  of  Equation  D15  gives 


2 C 


( w,  - wt ) 


(015) 


(016) 


1 

I 


and  a second  differentiation  gives 


JjT*  wi”a"  w*  ) 


(017) 
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Using  the  easily  derived  relations  that 

a*  + a • 2 £ (018) 

and 

a*  a"  « <u*  (019) 

the  solution  of  the  force  equation  can  be  easily  verified  by  substitution. 
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APPENDIX  E 

CONSTANTS,  PARAMETERS,  AND  CONVERSION  FACTORS 
1.  CONVERSION  FACTORS 


density 

1 kg/m3 

= 3.6127292  x 10~5  Ibm/in 

= 6.2427961  x 10~2  lbm/ ft 

energy 

1 joule 

= 0.737562  ft  lb 

length 

1 meter 

= 39.370039  In 

= 3.2808399  ft 

mass 

1 kg 

= 2.2046226  lbm 

1 grain 

= 6.479891  x 10'5  kg 

= 1.4285714  x 10~4  lbm 

pressure 

1 bar 

= 105  N/m2 

= 14.503774  lb/in2 

velocity 

1 km/sec 

= 1 m/msec  = 1 mm/psec 

= 39.370079  in/msec 


2.  CONSTANTS 

E = modulus  of  elasticity  in  tension.  Young's  modulus 
G = modulus  of  rigidity,  modulus  of  elasticity  in  shear 
v = Poisson's  ratio 


p = density 

E (kbar) 

2024  aluminum  sheet  724 

7075  aluminum  sheet  710 

JP-4 

JP-5 

water 


G (kbar)  ^(kg/m3)  C*  (m/s) 


276 

2770 

5112 

269 

2800 

5034 

785 

816 

1361 

1000 

1495 

* for  aluminum. 


C is  derived  from  C = 
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3.  CONSTANTS  FOR  SOME  TYPICAL  PROJECTILES 


Caliber  - Description 
(bullet  length  in  mm) 

Projectile 
Length  (mm) 

Diameter 

(mm) 

Mass 

(grams) 

Muzzle 
Velocity 
(km/  sec) 

7.62  - Ml 943,  Type  PS  (39) 

26.97 

7.87 

7.97  - 
8.035 

0.707  - 
0.732 

7.62  (54) 

29.97 

7.89 

9.59  - 
9.655 

0.8696 

30  cal  AP  (M2)  (63) 

35.51 

7.82 

10.74 

0.8382  - 
0.8534 

12.7  API  (108) 

63.5 

12.97 

48.84 

0.8382  - 
0.8687 

50  cal  AP  (M2)  (99) 

58.67 

12.97 

45.55 

0.8543 

14.5  Heavy  API  (113.8) 

50.55 

14.9 

64.22 

0.9997 

23  API  (150.5) 

106.5 

22.56 

187.8 

0.8839  - 
0.9144 

■i 
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